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FORTRAN PROGRAM FOR COMPUTING COORDINATES OF CIRCULAR ARC 
SINGLE AND TANDEM TURBOA/IACHINERY BLADE SECTIONS ON A PLANE 
by William D. McNally and James E. Crouse 
Lewis Research Center 

SUMMARY 

A FORTRAN IV computer program is presented which computes and plots coordi- 
nates for circular arc blade sections on a plane. Either single blade sections or tandem 
blades sections with up to 5 segments per blade section can be designed. Surfaces of 
blade segments consist of sir^le circular arcs. The arrangement of blade segments 
with respect to each other (for tandem blades) depends on the input parameters that 
specify gap, overlap, and convergence between the segments. 

Input is brief and can be altered rapidly. Input parameters describing the overall 
blade section include chord, camber, solidity, and inlet blade angle. Input to describe 
individual segments of the blade section include chord, camber, gap between adjacent 
segment and local segment, overlap of segments, maximum segment thickness, and 
radii of segment leading- and trailing-edge circles. Output consists of three main parts: 
(1) coordinates of individual segments suitable for making machine drawings, (2) geo- 
metrical input for companion blade-to-blade ideal flow programs, and (3) a Calcomp 
plot of the computed blade section in cascade at the input blade angle. All parts of the 
program except the plot routines are in general FORTRAN IV code and could be easily 
transferred to other IBM equipment. The plot routines, a short but important part of 
the output procedures, use a NASA Lewis code and would require recoding for use on 
other equipment. 

This report includes a listing of the FORTRAN IV computer program, with an ex- 
planation of the input required and the output generated. Numerical examples are also 
included. Running times are about 1/4 minute per data set on IBM 7094 equipment. The 
report does not include derivation or explanation of the equations on which the prc^ram 
is based. 



INTRODUCTION 


Specialized airfoil shapes are needed for todays highly loaded, high-speed compres- 
sors and turbines to avoid choking and premature separation. Shapes under study in- 
clude sii^le segment airfoils, airfoils with slots, and multiple segment airfoils in a 
tandem arrangement. 

Many of the single and tandem blade designs being studied have airfoil surfaces con- 
sisting of single circular arcs. The computation of geometry for such airfoils, particu- 
larly when placed in a tandem arrai^ement with controlled slot parameters, is compli- 
cated by the geometric calculations. 

This report describes a computer program for generating coordinates for circular 
arc airfoil shapes. One blade section is designed for each set of input data. A blade 
section consists of one cut through a blade at a given radius from the axis of rotation. 

The blade section may be composed of just one segment , or it may be a tandem section 
with two to five segments . The arrangement of blade segments with respect to each other 
(for tandem blades) depends on input parameters that specify gap, overlap, and conver- 
gence between the segments. The program does not provide radial stackii^ of blade 
sections, since only one section is designed for each set of data. 

Input is brief and can be prepared quickly. It consists entirely of geometric param- 
eters describii^ the overall blade section and the individual blade segments. Output 
consists principally of blade coordinates usable in other pregrams for the study of ideal 
flow and boundary layer. Output also includes coordinates usable for drafting or ma- 
chining, as well as a view of the blade in cascade in the form of a Calcomp plot. 

One of the principal uses of such a program is in conjunction with other computer 
programs for the analytical study of the performance and flow through turbomachine 
blading. This pregram permits the user to quickly generate and visualize circular arc 
blade shapes. The procedures of references 1 to 4 are then used to calculate velocities 
and streamlines on blade-to-blade stream surfaces of selected designs. The pregram of 
reference 5 calculates boundary-layer parameters from known flow velocity distribu- 
tions, and finally a program based on reference 6 calculates turbomachine losses from 
boundary- layer parameters at the blade trailing e(^e. These programs give the eigi- 
neer the ability to investigate blade shapes by testing them analsdically in a computer 
experiment. 

This report includes a listing of the program and a description of its input and out- 
put. The development of equations for the program is lengthy and will not be included. 
Internal pregram variables are not defined unless they are part of the input or output. 
Numerical examples are included to illustrate typical input values and the form in which 
output is given. 
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SYMBOLS 


C blade segment chord (fig. 3), ft; m 

F ratio of gap at inlet of channel between blade segments to gap at outlet of channel 
(figs. 3 and 12) 

G gap between blade segments (figs. 3 and 12), ft; m 
L overlap between blade segments (figs. 3 and 12), ft; m 
R radial coordinate direction (fig. 2) 

radius from axis of rotation to plane of blades (fig. 2), ft; m 
RI leading- edge radius of blade segment (fig. 3), ft; m 

RO trailing-edge radius of blade segment (fig. 3), ft; m 

S blade -to-b lade spacing on the cylindrical surface (figs. 1 and 2), ft; m 

TC total chord of overall blade section (fig. 1), ft; m 
TM maximum thickness of blade segment (fig. 3), ft; m 
Z axial coordinate direction (figs. 2 and 3) 

Ax overall blade section camber (fig. 3), deg 

inlet blade angle with respect to Z axis (fig. 3), deg 
d tangential coordinate direction (fig. 2) 

a solidity (fig. 1), TC/S 

<p camber of individual blade segment (fig. 3), deg 


GENERAL DESCRIPTION OF PROGRAM 
Characteristics of the Program 

From given inputs, the program calculates and plots coordinates of either a single 
blade section (see fig. 1(a)) or a tandem blade section (f^. 1(b)) with up to five segments 
per blade. (The plane of fig. 1 is the unwrapped cylindrical surface of fig. 2. ) 

All surfaces of the generated blade segments are single circular arcs tai^ent to 
leading- and trailing-edge circles. The radii of these arcs are a function of blade seg- 
ment cambers, chords, and thicknesses, which in turn are functions of given input pa- 
rameters. The position of the blade segments with respect to each other is also a func- 
tion of the inputs. 
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Figure 3. - Input variables. 


The input parameters (fig. 3) describe both the overall blade section and the individ- 
ual blade segments. The overall blade section is specified by a chord TC, camber Ak, 
solidity a = TC/S (see fig. 1), inlet blade angle and radius from axis of rotation to 
cylindrical surface (see fig. 2) . Individiaal blade segments are described by ratios 
of segment chord to the chord of the first blade segment C/C(l), ratios of segment cam- 
ber to first blade segment camber <p/<p(l), and ratios of maximum thickness and 
leading-edge and trailing-edge radii of the segment to local segment chord TM/C, RI/C, 
and RO/C. Segments are related to each other by the gap between them (in the form of 
ratio to total chord, G/TC), their overlap (also a ratio, L/TC), and the convergence in 
the channel between them F. (The chord, camber, gap, overlap, and convergence in- 
puts for the blade segments are not used when the blade section consists of only one seg- 
ment. ) 

For a tandem blade (more than one segment), the program follows an iterative pro- 
cedure in order to properly size the segments in relation to each other. From total 
camber Ak and total chord TC and the segment camber ratios cp/(p{l) and chord ra- 
tios C/C(l), initial estimates of segment cambers <p and chords C are calculated. 
Circular arc centerlines are fitted to these chords and cambers. The surfaces are also 
circular arcs that are tangent to leading- and trailing-edge circles, and meet the maxi- 
mum thickness requirement. Finally, the segments are located with respect to each 
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other. At this point the total camber formed from all estimated parameters is computed 
and checked against input total camber Ak. Adjustments are made, and the entire pro- 
cedure repeated until convergence is reached on total camber. After convergence, blade 
section coordinates and other output parameters are computed, and a plot of the blade 
section is made. 

Output from the program consists of printed computer listings and a Calcomp plot. 
The computer listings are divided into two main parts: (1) surface coordinates of indi- 
vidual blade segments suitable for making machine drawings and (2) geometrical input 
for blade-to-blade ideal flow programs (refs. 1 to 4) or a boundary layer program 
(ref. 5). The Calcomp plot shows the generated blade row at the input blade angle. Two 
overall blades are plotted with the proper solidity in order to identify the flow passage. 

The program is run at NASA Lewis on the IBM 7094-7044 direct coupled system with 
a 32 767 word core (77777^0^). The total program storage requirement is 65403^g^ of 
which 31717^gj is used in the storage of variables. The program runs in about 1/4 min- 
ute per data set on IBM 7094 equipment. 


Limitations of the Program 

The following are the principal limitations of the program: 

(1) Blade sections are generated on a plane surface, rather than a conical or merid- 
ional flow surface which would be more closely alined with the streamline flow when 
there is significant streamline slope. 

(2) Blade segment surfaces are single circular arcs. Multiple circular arcs or 
other types of variable geometry are not calculated by the program. 

(3) Each set of input data generates only one blade section. The program does not 
provide radial stacking of blade sections after several sections have been run. 

(4) The plotting portions of the program use routines that were developed at 
Lewis and would not be available or would need modification before they could be used 
on other machines. All other parts of the program, however, are in FORTRAN IV code, 
and could be easily transferred to other IBM equipment. 


Use of Program 

At Lewis, the program is being used to define blade sections for analytical para- 
metric studies using the programs of references 1 to 6. The Calcomp plots allow pre- 
liminary screening of cascades formed by applying the input variables over a wide range. 
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Selected configurations are then examined anal 5 dically for ideal flow, boundary-layer 
development, and losses. Some of these sections are later selected for experimental 
study. 

For applications in two-dimensional cascades or where radius does not change 
much across blade sections, output can be used for fabrication purposes. 


NUMERICAL EXAMPLES 

Two numerical examples are given which illustrate the use of the program. The 
first is a two- segment tandem blade section, and the second is a three-segment blade 
section with the front section acting as a slat. Both blade sections are designed for the 
same overall parameters which are listed in table I. The input for these two examples 


TABLE I. - OVERALL DESIGN PARAMETERS 
FOR TWO- AND THREE-SEGMENT 


TANDEM BLADE SECTIONS 


Total chord, TC, ft 

0. 18583 

Solidity, a 

1. 235 

Overall camber, Ak, deg 

72.24 

Inlet blade angle, deg 

56. 53 

Radius from axis or rotation, R^, ft 

0. 77080 


and the generated plots of blade shapes appear in figures 4 and 5. These examples il- 
lustrate typical values of input parameters for two- and three-segment blade sections. 
Sample output for the first example is discussed under OUTPUT. 


INPUT 

Figure 6 shows the placement of input variables on data cards. The first input card 
is for a title which identifies the data set and is printed on the output. The user may 
type whatever information he wishes in any of the first 72 columns of this card. The 
remaining cards are for input data. The input variables are defined in the next section. 
Further explanation of the proper preparation of input is contained in the section Typical 
Values and Limits of Input Variables. 
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Circumferential distance Circumferential distance 


Axial distance 



EXAMPLE - TwC SBCTICN TANChM BLADE 


N TChCRO SCUD CELK 

2 0.iB563 1.235C 72,2A0 

C/C(l) ARRAY 

1 -25CCO 

PHI/PHI(1> ARRAY 




1 .56250 

G/TC 

ARRAY 

0.025CO 

L/7C 

ARRAY 

0. lOOOO 

f ARRAY 

1 .^CCCO 

7M/C 

ARRAY 
0. 120CC 

0. ICCCO 

RI/C 

ARRAY 
O.C 15CC 

0.C2CCC 

RO/C 

ARRAY 

0.00800 

C. 0C7OO 


Figure 4. - Input and generated plot of two-section tandem blade example. 


KAPIN 

56.530 


example - THREE SECTiC^ 7A^CEM BLADE 

N rCHCkO SGLIO CELK KAFIK 

3 0.165E3 1.235C 72.2A0 56.530 


C/U 1 ) ARRAY 



5 .ICCOO 

6. CCCCO 

FHI/PHin) ARRAY 

0.5C0C0 

C. 60CCC 

G/TC ARRAY 

C . C125C 

C. 02ACC 

L/7C ARRAV 

0.CA2CC 

0. 057CC 

F ARRAY 

1.2CCC0 

1 . ACCCO 

7M/C ARRAY 

0. J4C00 

0. 12C00 

C. 12CCC 


RI/C ARRAY 

O.OeiOO O.C3200 C.C25CC 


RO/C ARRAY 

O.OiOCC 0.CC700 0.CC700 

Axial distance 


Figure 5. - Input and generated plot of three-section tandem blade example. 




RAUUS 

0.77080 


R ACIUS 
0-770B0 
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Figure 6. - Input data form. 


Input Variables 

Schematic representations of these variables appear in figures 1 to 3. After the 

title card, the following input variables are given: 

N integer number (1 to 5) of blade segments comprising the blade section; 

equals 1 when designing a single, circular-arc blade section; must occupy 
column 10 of the data card (fig. 6) 

TCHORD total chord of the overall N- segment tandem blade section TC, ft; m 

SOLID solidity of the blade row, a, that is, total chord divided by blade spacing 

TC/S. (Solidity is only used in the plotting part of the program to produce 
a duplicate blade on the plot. ) 

DELK total camber of the overall blade sections, A/c, deg 

KAPIN blade inlet angle or angle between tangent to mean camber line at leading 
edge of first blade segment and the Z axis, deg 

RADIUS radius from axis of rotation to cylindrical blade plane R^, ft; m (RADIUS 
is only used to convert tangential coordinates, R0, in feet or meters, to 
radians for input to the ideal flow programs, refs. 1 to 4.) 
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Each of the following arrays has N - 1 entries. K N = 1, a blank card should be 
given for each of these 5 arrays. 


COCl array of ratios of chords of blade segments 2 to N to the chord of the first 

segment, C/C(l) 

PHOPHl array of ratios of cambers of blade segments 2 to N to the camber of the 
first segment, <p/cp{\) 

GOTC array of ratios of gaps between blade segments to the total chord of the over- 
all blade section, G/TC 


LOTC array of ratios of overlap between blade segments to the total chord of the 
overall blade section, L/TC 

F array of channel convergences between blade segments, F (F is the ratio of 

the gap at the channel inlet to the gap at the channel outlet. ) 

Each of the arrays below has N entries, one for each of the blade segments; 


TMOC array of ratios of maximum blade segment thickness to chord of the individual 
blade segments, TM/C 

RIOC array of ratios of leading- edge radius to chord of the individual blade segments, 
RI/C 

ROOC array of ratios of trailing-ec^e radius to chord of the individual blade segments, 
RO/C 


Typical Values and Limits of Input Variables 

Ranges of typical values are given in this section for the input variables. Limits 
are also given beyond which unreasonable blade sections (and hence errors in the pro- 
gram) will occur. 

N, the number of blade sections, can be any integer from 1 to 5. For typical tandem 
blades, N is usually 2 or 3. To design a single blade section, N is set equal to 1, 
and blank cards are used for the COCl, PHOPHl, GOTC, LOTC, and F arrays 
(Fig. 7 is the input and the corresponding output plot of a single blade section. ) 

Since N is an integer, it must be right shifted on the data card; that is, it must 
occupy column 10 (see fig. 5). 

TCHORD can be any positive value. 

SOLID can also be any positive value; the range from 0. 5 to 2. 0 is typical. 
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Circumferential distance 


Axial distance 


sample blade wnh a single segment 

N TCHCRD SOLID CELK KAPIN RADIUS 

1 0.1S20C 1.3CCC 50-000 ^5.000 0.62500 

C/C( 1 ) ARRAY 

- 0 . 

FHI/PHim ARRAY 
- 0 . 

6/TC ARRAY 

- 0 . 

L/TC ARR^y 

- 0 . 

F ARRAY 

- 0 . 

TM/C ARRAY 

0. 11000 

RI/C ARRAY 

0.01300 

RfJ/0 ARRAY 

0 . 00 8 0 0 

Figure 7. - Input and generated plot of blade section with single segment. 

DELK, the overall chamber, must be a positive number or zero. Values as high 
as the 180° will run, but the range from 5° to 120° is typical. If DELK has a small 
value (from 0° to 10°) the program will not converge to an answer if other param- 
eters such as segment camber, gap, overlap, and convergence are not physically 
compatible with DELK. 

ICAPIN, the blade inlet angle, can be positive, negative, or zero. Values between 
the limits of -90° to 90° are allowable, but the range from -30° to 70° is typical. 

RADIUS can be any positive value. 

TCHORD and RADIUS are the only inputs with units of length; units should be the 
same on these two variables. Generally either feet or meters are used so that output 
can be used with the ideal flow programs (refs. 1 to 4). This is not required, however, 
and any units of length are acceptable. Units on all output coordinates will always cor- 
respond to what was used on these two input quantities. 

COCl can be any positive value. The range from 0. 1 to 10. 0 is typical. 
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PHOPHl can be any positive or negative value, or zero. Values from 0 to 3. 0 are 
most common. To obtain a very straight front segment, PHOPHl should contain 
very large values. To obtain a very straight aft segment, PHOPHl should be near 
or equal to zero for that segment. 

GOTO can be any positive value from zero to about 0. 5 depending on other inputs 
such as segment cambers, overlaps, and convergences. The range from 0. 01 
to 0. 04 is typical. 

LOTC can have positive or negative values, or be zero. Typical values are con- 
tained in the range from 0. 1 to -0.05. Values above 0.4 or below -0.2 will gener- 
ally cause errors and prevent the program from running, 

F can have positive values from zero to about 10. 0. The range from 0. 9 (diverging 
passage) to 1. 5 (converging passage) is most typical. When F = 1,0, the capture 
area of the passage between blade segments is equal to the exit area of the passage. 

TMOC is allowed positive values from zero to about 0, 8. Values in the range from 
0. 1 to 0. 2 are most typical. Elements of TMOC must be at all times at least twice 
as large as the corresponding elements of RIOC and ROOC in order for the program 
to run. (K TMOC equals zero, RIOC and ROOC must also be zero for that blade 
segment. ) 

RIOC and ROOC may have positive values from zero to about 0. 4. Most values are 
in the range 0. 01 to 0.1, Corresponding elements of RIOC and ROOC do not have 
to equal each other. (RIOC may only equal zero if the corresponding element of 
ROOC also equals zero. ROOC, on the other hand, may equal zero at any time, 
regardless of the values in RIOC. ) 

Example of Adjustment of Inputs In Design Process 

The program is used here to design a two- segment tandem blade section. Given 
the overall blade section parameters, an initial selection is made for the other input 
variables. These variables are subsequently changed (twice in this example) until a 
final blade section is accepted. 

Changes are made after inspection of the machine plots which accompany the com- 
puter output. They are made to obtain a blade section which appears to have a good flow 
path while satisfying the overall blade parameters. These iterations on input variables 
also illustrate the effect of the different input parameters on the final blade shape. 

The blade section to be designed has the overall blade parameters listed in table II. 
In order to obtain an initial picture of a blade section meeting these specifications, gen- 
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TABLE n. - OVERALL DESIGN 


PARAMETERS FOR TWO- 
SEGMENT TANDEM 


BLADE SECTION 


Total chord, TC, ft 

0. 192 

Solidity, a 

1.3 

Overall camber, A/c, deg 

50.0 

Inlet blade angle, deg 

45.0 

Radius of rotation, R^, ft 

0.625 


TABLE in. - VARIABLE INPUT PARAMETERS FOR TWO-SEGMENT TANDEM BLADE SECTION 


Run 

Ratio of seg** 

Ratio of seg- 

Ratio of 

Ratio of 

Ratio of gap 

Ratio of maxi- 

Ratio of leading- 

Ratio of trailing- 


ment chord to 

ment chamber 

gap to 

overlap 

at channel 

mum thickness 

edge radius to 

edge radius to 


chord of first 

to first blade 

total 

to total 

inlet to gap 

to local seg- 

local segment 

local segment 


blade segment, 

segment 

chord, 

chord, 

at channel 

ment chord, 

chord, 

chord. 


C/C(l) 

chamber, 

G/TC 

L/TC 

outlet, 

TM/C 

RI/C 

RO/C 



cp/(p{l) 



F 




1 

1.0 

1.0 

0.05 

0. 10 

1. 5 

0. 15 

0. 02 

0.01 







. 15 

. 02 

. 01 

2 

1. 3 

1.6 

0. 03 



1. 1 

0. 13 



0.007 







. 13 


.007 

3 









0. 11 

— 





. i 


1. 2 j 

. 12 

— 



eral initial values of the other input parameters were chosen and run. These values are 
listed in table III (run 1). The resulting blade section is shown in figure 8(a). 

Changes made after an initial run on the program are entirely based on the user’s 
experience and his concept of the final desired blade shape. For this example we wanted 
more chord and camber to be concentrated in the rear blade segment; so, in run 2, 
C/C(l) was increased from 1.0 to 1.3, and cp/ip (1) from 1.0 to 1.6. The channel gap 
was also decreased (G/TC = 0. 05 to 0. 03), as well as the channel convergence between 
blade segments (F = 1. 5 to 1. 1) in order to bring the segments closer together. Finally 
the blade thicknesses TM/C and the outlet radii RO/C were reduced. The blade sec- 
tion resulting from run 2 (table III) is pictured in figure 8(b). From experience it ap- 
peared that this blade section was still thicker than desired and that its channel needed 
more convergence. Appropriate changes were made for rim 3 (table III), and the final 
blade section is shown in figure 8(c). This section was accepted for further analysis by 
the ideal flow programs (refs. 1 to 4). 
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OUTPUT 


Output from the program consists of two principal parts: a computer listing with 
printed tables of output variables and a Calcomp plot that pictures schematically the 
generated blade. 

A sample computer listing for the two- section tandem blade example is given in 
table IV. In this table some sections of the output have been abbreviated because they 
were too long. In all cases output labels agree with program variable names which are 
defined in the next section. 


TABLE IV. - SAMPLE OUTPUT FOR TWO-SECTION TANDEM BLADE EXAMPLE 


example - TWO SECTION TANOEM BLADE 


N TCHORD SOLID OELX KA F I N RADIUS 

2 0.16583 1.235C 72.240 56.530 0.77080 

C/C( 1 » ARRAY 

1.25CCC 

FHI/PHUI) ARRAY 

1.56250 

&/TC ARRAY 

C.C25CC 

L/TC ARRAY 

t . ICCCO 

f ARRAY 

1 . 4CCCC 

IM/L ARRAY 

C. 12000 O.IOOOO 

Rl/L ARRAY 

C.015CC 0. 02000 

RO/C ARRAY 

c.ooaoo 0.00700 


LVtRALL liLACE PARARtltRS 

N KHORO PITCH bCLlC CELK KAPIN GAM8 THUS XD8 YCb 

2 C. 1 8563 C. 1 5047 1.2350 72.240 66.5<00 19.6156 C.162^^ 0.17413 


BlALE itcMEM NO. I 


CHURO 

RI 

RO 

IHETA 



C .CS 334 

C. LC14C 

C. CCC75 

C. 41051 



X i 

Y1 

XQ 

YC 

XCM 

Y CR 

C .C014C 

C.C014C 

C. C5259 

C. 000 75 

0.04534 

0 .0 0 80 3 

XI 

Y 1 

X2 

Y2 

GAR 

GAf^R 

C . 

-c. 

0. C9 334 

-C. 

0 . 

0. 

Fhl S 

RS 

HS 

BS 

KIS 

KOS 

5 6 .52 24? 

C.CS737 

C. C457C 

C.0837 3 

2 7 .48966 

29.0 3 381 

FHIC 

RC 

hC 

ec 

K IL 

KOC 

34 .70661 

C. 1528 ? 

C. C4595 

C. 14483 

16.94380 

17.76481 

PHIP 

RP 

hP 

BP 

KIP 

KOP 

12 .24C54 

C. 42655 

C. C**7C1 

C. 42416 

6. 1 1793 

6.12261 

XX 

YS 

YP 

NCtL = 47 



C . 

C. CC225 

-C. CCC16 




C .C02CC 

C. CC326 

C, CCC06 




C .CC4CC 

C.CC426 

C.CCC26 




C .C06C0 

C.CC518 

C. CCC46 




C .CC6CC 

C.CC6C4 

C. CCC65 




C .CICLC 

C.CC686 

C. CCC83 




C.C12CC 

C.CC762 

C. CClOO 




C .Cl <CC 

C. CC633 

0.00116 
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TABLE IV. - Continued. SAMPLE OUTPUT FOR TWO-SECTION TANDEM BLADE EXAMPLE 


0.01600 

C.C0900 

0.00131 

O.CIEOO 

0.00961 

0.00145 

C.02000 

C.01018 

0.00158 

C.C2200 

C. 01071 

O.CC170 

C.024G0 

C.01119 

0.00182 

C.0260C 

€.01162 

0.00192 

€.02600 

C.C1201 

O.C0201 

C. 03000 

0.01236 

0.00210 

C. 03200 

€.01267 

0.00217 

C .03400 

C. 01293 

0.C0224 

€.03600 

C.01315 

0.00230 

0.03600 

€•€1333 

0.00234 

C. 04000 

€.01347 

0.00238 

C.C4200 

0.01357 

0.00241 

€.04400 

€.01362 

0. 00243 

C.0460C 

€.01364 

0. 00244 

C.C460Q 

€.01361 

0.00244 

C.05CCO 

C.C1354 

0.00243 

C.C52CQ 

€^01343 

0.QQ241 

C.C5400 

C. 01328 

0.CC23 8 

C.0S600 

0.01309 

0.00234 

C.C580C 

C.CI286 

O.C0230 

C.C6000 

C.C1258 

0.00224 

C.C62C0 

0.01226 

0.0C217 

0.06400 

C.C1190 

0.00210 

C.C660C 

€.€1150 

0.00201 

C.C660C 

C. 01105 

C. 00192 

C.C7CC0 

0.01055 

0.00182 

C.C72GC 

C.010C2 

0.CC171 

O.C7400 

0.00943 

0.CC158 

C.C76CC 

C.C0880 

0.00145 

C.C7600 

C.CC8I2 

0.00131 

C.C8CC0 

C.CC739 

0.00116 

0.08200 

0.00662 

0.00100 

C.C8400 

C.CC579 

0.00083 

C.C86CC 

C.CC49C 

€.00065 

C.C8600 

C.C0397 

0.0CC46 

C.C9C0C 

C.C0297 

0.CQC27 

0.C92Q0 

€.00192 

0.0CC06 


8LACE SE6HEKT NO. 2 


CI-43RD 

R1 

RO 

THETA 



C. 11667 

C.0C233 

0.€0C82 

0.76547 



XI 

71 

XO 

YC 

XCN 

YCM 

O.C0233 

C.C0233 

0.11586 

0.00082 

0.05386 

0.01522 

XI 

Y1 

X2 

Y2 

GAH 

GAMR 

C.C7C86 

C. 00825 

C. 17448 

0.06189 

27.36919 

27.36919 

PHIS 

RS 

HS 

BS 

KIS 

KOS 

7C.04290 

C. 10049 

0.C5669 

0.07939 

33.62912 

36.41378 

PHIC 

RC 

HC 

BC 

KIC 

KOC 

54.23220 

C. 12454 

0.05761 

0.10927 

26.35063 

27.88157 

FHIP 

RP 

HP 

BP 

KIP 

KOP 

37.66516 

0.17337 

0. €5922 

0.16391 

18.89117 

18.97399 

G 

6A 

GAOC 

L 

F 

FA 

€.00465 

C. 00463 

0. 04960 

0.01858 

1.40000 

1.40619 

XX 

YS 

YP 

NOEL => 24 



0. 

C. 00357 

-0.00096 




C.005C0 

C.C0676 

0.C0C77 




C.OIOOC 

C.C0959 

C. 00233 




C .01500 

0.01204 

0. 003 73 




C.C20CO 

0.01416 

0.00497 




C.C25C0 

0.01597 

0.00605 




C.C3C00 

€.01748 

0.00698 




C.C3500 

C. 01672 

0.00776 




COMPITED ! 

INPUT FOR IDEAL FLOW 

PROGRAMS 




BLAD€ 

MCHORO 

STGR 

RSTGR 

R1 

RO 

1 

0. 07284 

0.C7474 

0.05761 

0.00140 

0.00075 

2 

O; 11442 

0.02924 

0.02254 

0.00233 

0.00082 


BLADE 

BETIS 

BETOS 

BETIP 

BETCP 

1 

67.07586 

10.55239 

45.70413 

33.46359 

2 

45.84613 

-24.19677 

31.10318 

-6.75698 


SINC 

-9.55789 


ML 6 

0. 

0.06054 

THLE 

o.ocooo 

0.C5I37 

RTHLE 

0.00000 

0.03960 

MTE 

0.07284 

0.17496 

THTE 
0. 07474 
0.08062 

RTHTF 

0.05761 

0.06214 

MCL 

0.00140 

0.06287 

THCL 

O.GOOOC 

0.05137 

RTHCL 

c.ooooo 

0.03960 

MCT 

0.07209 

0.17415 

THCT 

0.07474 

0.08062 

RTHCT 

0.05761 

0,06214 
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TABLE IV. - Concluded. SAMPLE OUTPUT FOR TWO- SECTION TANDEM BLADE EXAMPLE 


BLAtE 5EOMEM NQ. i 


MSPS 

IHSPS 

BTHSPS 

-{.CCC22 

-C.00C31 

-C,CCC24 

( .tCC66 

O.OC237 

C. 00183 

(.CC15E 

C. CC5CC 

0.CC386 

C.OC254 

C. CC758 

0. 00584 

C.OCi53 

C.OICIC 

0, CC778 

C.CC455 

C. 01256 

C. CC966 

( .OC561 

0.01496 

C. 01155 

( .CCt69 

0,01735 

C. 0133 7 

c.oc?ai 

C. C1966 

C. 01516 

C.CCfc96 

0.02193 

0.01691 

C.O iOlA 

C.C2416 

0.01862 

C.O 1 

C. 02634 

C. C203C 

C.O I25ti 

0.02847 

C. 02194 

C.0138A 

C, C3056 

C. 02355 

C .Olf 14 

0. C326C 

0. C2513 

C.01t4£ 

0.C346C 

C. 02667 

(.01 ?ec 

C. 03656 

0. 02818 

C.G ISU 

0.03848 

0. 02966 

c.oicss 

C.04C35 

C. 0311C 

C.022CC 

C, C4219 

0. 03252 

C ,02246 

0. 04398 

C. 0339C 

C.0<494 

0.C45 73 

C. 03525 

C .0 2 t44 

C. C4744 

C. 03656 

C .02 798 

0. 0491 C 

C. 03785 

C .02953 

C.C5073 

0. (391 C 

C.O 2 112 

0.05232 

C. C4033 

C .0 227 3 

C,C5386 

0. 04152 

( .0243 7 

C. 05536 

0. 04266 

C. 02(03 

0.C5663 

0,04380 

C .C 2 7 72 

C. 0582 5 

0. C449C 

C - 0 2 9 4 4 

C. C5962 

C. C4596 

(.04118 

C. 06096 

U. 04699 

( .04295 

C.C6225 

C. 04798 

C. 044 75 

C. C635C 

C, 04895 

C. 04(56 

C. 06471 

0.04987 

( .0^643 

C. 0656 7 

0.05077 

( .050 32 

C. C6698 

0. 05163 

(.C5^^ ? 

C. 06805 

0. 05245 

C .05416 

C. C690 7 

0.05324 

(.05(15 

C.C70C5 

0.05399 

(.05(16 

C. 0 709 7 

C. 054 71 

( .C6C 19 

C. 071 85 

0. 05536 

( .0 ( 2 2 6 

3. 0 726? 

0. 056C2 

( .( (^2 7 

C.C 7 344 

C. 05661 

C ,C (( 5C 

C.C7416 

0.05716 

( ,((666 

0.07482 

C. 05767 

( .C 7(89 

C. C7542 

0.05814 


6LACE EtOt'tM 

NQ. 2 



IHSPS 

BTHSPS 

-( .0 CC21 

C.CCG93 

0. 00072 

( . 30 4JC 

C-0C637 

0.00491 

0 .UC629 

C. Oil 30 

C. 008 7 1 

( .0 i^et 

0. 015 7e 

C. 0121 6 

( . 0 1 71 0 

C. C1984 

0.01529 

( ,C21(C 

0.02351 

0,01812 

(.Ji(l 7 

0. 02681 

O.C2066 


PSPP 

ThSPP 

RThSPP 

0.00 132 

-0.00272 

-0.CO21C 

0.00272 

-0.00085 

-0.CC065 

0.00413 

0.00101 

0.00078 

0,00554 

0.00286 

0.00221 

0.00696 

0.00471 

0.003(3 

0.00039 

0,00654 

0.0C5C4 

0,00983 

0.00836 

0.00(44 

0,01126 

0.01017 

0.00784 

0,01271 

0.C1198 

0.00923 

0.01416 

0,01377 

C.C1062 

0.01562 

0.01556 

0.C1199 

0,01708 

0.01733 

0,C133( 

0.01855 

0.01910 

0.01472 

0.02003 

0 .02086 

0.C16C8 

0,02151 

0,02260 

0.01742 

0.02300 

0.02434 

0,C187( 

0.02449 

0.02607 

0.C20 IC 

0.02599 

0.02779 

0.C2142 

0.02749 

0.C2950 

0.02? 74 

0.02901 

0.03120 

0.C24C5 

0-03052 

0.03289 

0.02535 

0.03205 

0.03457 

0.02(65 

3.03358 

0.C3625 

0.02794 

0 .0351 1 

0.03791 

0.C2922 

0.03665 

0.03956 

0.C3C49 

0.03820 

0.04121 

0.C31 7( 

0.03975 

0.04284 

0 . 0 3 3 C 2 

0.04 131 

0.04447 

0.03427 

0.04288 

0.046U8 

C. 03552 

0 .04445 

0.04769 

0.036 76 

0.04603 

0.04929 

0.03799 

0.04761 

0.05087 

0.03921 

0.04920 

0.05245 

C .04043 

0,050 79 

0.05402 

0.041(4 

0.05239 

0.05558 

0.04284 

0.05400 

0.05713 

0.C44C4 

0.05561 

0.05867 

0.C4522 

0 .05723 

0. 06020 

0.04(40 

0 .05886 

0.06173 

0,047 5 6 

0.06049 

0.06324 

0. 04 8 74 

0.062 13 

0.06474 

0.04990 

0 .06 37 7 

0.06623 

'1.C51C5 

0 .06542 

0.06772 

O.C522C 

0.06707 

0.06919 

0.05333 

0.06874 

0,07066 

0 .0544( 

0.07040 

0.07211 

0. 05559 

0.07208 

0.07356 

0.C567C 


PSPP 

ThSPP 

RTHSPP 

0 .00075 

-0.00482 

-0.00372 

0.00527 

-0.00125 

-C.0CC97 

0.00983 

n. 00210 

0.CC162 

0.01442 

0.00525 

0.CC4C4 

0.U1904 

0.00819 

0.006 3 1 

0.02370 

0.01094 

O.CC843 

0,07839 

0.01349 

0 . ac4c 


Output Variables 


The first page of output contains a copy of the input to the program. These variables 
are defined in the Input Variables section on page 9. 

The next page of output lists some overall blade section parameters, some of which 
are repetitions from the input list. The others are defined as follows: 

PITCH blade-to-blade spacing S in the 6 direction or ratio of total chord to 

solidity TC/a (fig. 1), ft; m 

GAME angle between chord line of first blade segment and chord line of overall 

blade section (fig. 9), deg 
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Figure 9. - Output variables for overall blade section. 


THETB angle between chord line of overall blade section and a line joining leading- 
edge circle center of first blade segment and trailing- edge circle center 
of final blade segment (see fig. 9); positive if RI(1) > RO(N) and nega- 
tive if RI(1) < RO(N), deg 

XOB(YOB) distance from first segment leading-edge line (first segment chord line) to 
circle center at trailing edge of final blade segment (fig. 9), ft; m 

Following the overall blade parameters are lists of parameters for each of the individ- 
ual blade segments: 

CHORD chord length of blade segment, that is, distance from leading-edge 

point to trailing-edge point (fig. 10), ft; m 

RI(RO) leading- (trailing-) edge radius of blade segment (fig. 10), ft; m 

THETA angle between chord line of blade segment and a line joining leading- 

and trailing-edge circle centers (fig. 10); positive if RI > RO, 
and negative if RI < RO, deg 

XI(YI) distance between leading-edge line (chord line) of blade segment and 

center of leading-edge circle (fig. 10), ft; m 
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Figure 10. - Output variables for individual blade segment 


XO(YO) 

XCM(YCM) 

Xl(Yl) 

X2(Y2) 

GAM 

GAME 


distance between leading-edge line (chord line) of blade segment and 
center of trailing-edge circle (see fig. 10), ft; m 

distance between leading- edge line (chord line) of blade segment 
and point on mean camber line at which slopes of both blade sur- 
faces are equal to slope of mean camber line (fig. 10), ft; m 

distance between leading- edge line (chord line) of first segment 
and leading- edge point of local segment (fig. 9), ft; m 

distance between leading- edge line (chord line) of first segment and 
trailing-edge point of local segment (fig. 9), ft; m 

angle between chord line of local blade segment and chord line of 
previous blade segment (fig. 9), deg 

angle between chord line of local blade segment and chord line of 
first blade segment (fig. 9), deg 


PHIC(PHIS, PHIP) overall camber of local blade segment mean camber line (suction 

surface, pressure surface) from a line through leading- edge 
circle center to a line through trailing-ec%e circle center (fig. 
11(a)), deg 

RC(RS, RP) radius of curvature of local blade segment mean camber line (suc- 

tion surface, pressure surface) (fig. 11(a)), ft; m 


HC(HS, HP) distance from leading- edge line of blade segment to center of cur- 

vature of mean camber line (suction surface, pressure surface) 
of blade segment (fig. 11(a)), ft; m 

BC(BS, BP) distance from chord line of blade segment to center of curvature 

of mean camber line (suction surface, pressure surface) of blade 
segment (fig. 11(a)); negative when PHIC( PHIS, PHIP) is negative 
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Figure 11. - Continuation of output variables for individual blade segment. 


When PHOPHl, and thus PHIC, of a segment equals 0, RC and BC are set to 
999.99999. 

KIC(KOC) angle between chord line of blade segment and tangent to mean camber line 
at the center of leading-edge (trailing- edge) circle (see fig. 11), deg 

KIS(KOS) angle between chord line of blade segment and tangent to suction surface at 
leading-edge (trailing-e<^e) transition point (see fig. 11), deg 

KIP(KOP) angle between chord line of blade segment and tangent to pressure surface at 
leading-edge (trailing- edge) transition point (see fig. 11), deg 

KIC(KIS, KIP) and KOC(KOS, KOP) are defined positive as shown in figure 11 for a 
centerline or surface which has positive camber. They will be negative for a surface 
with negative camber. 
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G 

GA 

GAOC 

L 

F 

FA 


gap between trailing ec^e of previous blade segment and suction surface of 
local blade segment (fig. 12); measured perpendicular to the chord line of 
previous blade segment along line passing through trailing- edge circle center 
of the previous blade segment, ft; m 

actual gap between trailing edge of previous blade segment and suction surface 
of local blade segment (fig. 12); measured perpendicular to suction surface 
of local blade segment along line passing through trailing-edge circle center 
of previous blade segment, ft; m 

ratio of GA to CHORD of previous blade segment (fig. 12) 

distance between gap G at trailing edge of previous blade segment and gap (FXG) 
at leading ec^e of local blade segment (fig. 12), ft; m 

ratio of gap FxG at leading edge of local blade segment to gap G at trailing e<%e 
of previous blade segment (fig. 12) FxG is measured perpendicular to chord 
line of previous blade segment along a line passing through leading e(^e 
circle center of local blade segment 

ratio of actual gap FAxGA at leading ec^e of local blade segment to actual gap 
GA at trailing edge of previous blade segment (fig. 12) (Actual gap FAxGA 
is measured perpendicular to a line (A -A in fig. 12) which bisects the tan- 
gents to the suction surface of the local blade segment and the pressure sur- 
face of the previous blade segment where the line FAxGA meets these sur- 
faces. The line containing FAxGA passes through the leading- edge circle 
center of the local blade segment. ) 



Figure 12. - Input and output variables in overlap region. 
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SINC angle between tangents to mean camber lines of local blade segment and pre- 
vious blade segment at points of intersection with line containing FXG (see 
fig. 12), deg (SINC is a measure of the incidence of the average blade-to- 
blade flow on the leading edge of the local blade segment. SINC is negative 
as shown in figure 12 since the mean flow would have negative incidence in 
this blade orientation. ) 

XX array of distances (parallel to blade segment chord line) between leading-edge 

line of blade segment and points at which blade surface coordinates (YS and 
YP) are given (fig. 10), ft; m 

YS(YP) array of perpendicular distances from chord line of blade segment to points 
on suction (pressure) surface of the segment (fig. 10), ft; m 

NDEL number of blade coordinate points (XX and YS, XX and YP) along suction or 
pressure surfaces of local blade segment. 

For blades with normal levels of positive camber, some values of YP at inlet and outlet 
may be negative. These are points on the pressure or suction surface circular arcs 
that occur prior to the leading-edge radius or after the trailing-edge radius (fig. 13). 



For a blade with small positive camber, or with negative camber, many values of YP 
(and sometimes YS) can be negative (fig. 14). 

Following the blade coordinates for each of the blade segments are output param- 
eters that serve as inputs for the ideal flow pregrams. These programs are reported 
in references 1 to 4. They compute ideal flow on an axisymmetric blade-to-blade sur- 
face of a single or tandem bladed turbomachine in either subsonic or mildly transonic 
flow. 

To obtain input to be used in the ideal flow programs, the flat plane in which the 
blade section lies is assumed to be wrapped about a cylinder of radius equal to the input 
parameter, (RADIUS, R^^, in fig. 2). This cylinder serves as the axisymmetric blade- 
to-blade surface required for the input of geometry to the ideal flow programs. 
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(a) Small positive or negative camber. 



(b) Negative camber. 


Figure 14. - Example blade sections with negative surface coordinate points. 

Specific output quantities which are required as geometric input parameters in the 

ideal flow prc^rams are defined in the following: 

MCHORD chord lengths of blade segments in Z direction (figs. 15 and 16), ft; m 

STGR angular 9 coordinates of trailing edges of blade segments with re- 

spect to leading edges of blade segments (figs. 15 and 16), rad 

RSTGR angular distances of trailing et^es of blade segments from leading 

edges of blade segments (figs. 15 and 16), ft; m 

RI(RO) leading-edge (trailing- edge) radii of the blade segments (figs. 10 

and 16), ft; m 

MLE(MTE) distances in Z- direction from leading edge of first blade segment to 

leading (trailing) edges of other blade segments (fig. 15), ft; m 

THLE(THTE) angular 9 coordinates of leading (trailing) edges of blade segments 

with respect to leading edge of first blade segment (fig. 15), rad 

RTHLE(RTHTE) ai^ular distances of leading (trailing) edges of blade segments with 

respect to leading edge of first blade segment (fig. 15), ft; m 

BETIS(BETOS) armies with respect to Z-direction at tangent points of leading- 

(trailii^) edge radii with suction surfaces of blade segments 
(fig. 16), deg 
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Figure 16. - Blade segment output variables used for plots and ideal flow programs (refs. 1 to 4). 



BETIP(BETOP) 


MCL(MCT) 


THCL(THCT) 

RTHCL(RTHCT) 


angles with respect to Z-direction at tar^ent points of leading- (trail- 
ing) edge radii with pressure surfaces of blade segments (fig. 16), 
deg 

distance in Z-direction from leading edge of first blade segment to 
centers of leading-edge (trailing- edge) circles of blade segments 
(fig. 15), ft; m 

angular 9 coordinates of centers of leading- edge (trailing- edge) cir- 
cles with respect to leading edge of first blade segment (fig. 15), rad 

angular distances of centers of leading- edge (trailing- edge) circles 
with respect to leading edge of first blade segment (fig. 15), ft; m 


The preceding variables are followed by the coordinates of suction and pressure sur- 
faces of the individual blade segments. These coordinates are given with respect to axes 
in the Z-direction passing through the leading-edge circle centers of each of the seg- 
ments. 


MSPS(MSPP) 


THSPS(THSPP) 


RTHSPS(RTHSPP) 


array of distances in Z-direction between leading edges of individ- 
ual blade segments and points on the suction (pressure) surface 
at which blade coordinates (THSPS and THSPP) are given as out- 
put (fig. 16), ft; m 

array of angular 9 coordinates from a line in the Z direction 
through the leading edge circle center of each segment, to points 
on the suction (pressme) surface of the segment (fig. 16), rad 

array of distances (RADIUS x THSPS(THSPP)) corresponding to 
THSPS(THSPP) (fig. 16), ft; m 


Output Blade Plots 

The Calcomp plot portion of the output shows the blade as it would appear in cascade 
at a given solidity and inlet blade angle. The plot is very helpful in evaluating the blade 
visually; the user can see immediately whether the shape resulting from the program 
agrees with his concept. 

The complete Calcomp plot for the two- section tandem blade example is shown in 
figure 17. All plots have the same format as this one. To the left of the plot are printed 
the complete input and some selected output variables. The input variables on the Cal- 
comp plot and in the program are related as follows: C/C(l) = COCl, PHI/PHI(1) = 
PHOPHl, etc. For the output, PHI are the blade segment cambers, PHIC. The blades 
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CHORD 
. ig'303 

COMBER 

72.2400 

soLion 

1.2360 

Y 

PITCH 

.ISOiiS 

KflPIN 

66.6300 

RHDIUS 

.77060 


C/Cil) 

1.2S00 

PHI/PHI ;i; 
1.6626 

) G/TC 
.0260 

L/TC F 

.1000 I. 4000 

IK/C 
.1200 
. 1000 

RI/C 

.0160 

.0200 

RO/C 

.0060 

.0070 


c 

.OG.331 

.11667 

pi-|] 

34 ! 7066 
64 2322 

GflM 

0.0000 

27.3691 

GflMR 
0.0000 
27. 3691 



"b.cc 


. )3 


Figure 17. - Full Calcomp plot for two-section tandem blade example. 



are drawn in a position corresponding to the blade angle, KAPIN. The Z-axis is nor- 
mal to the sides of the Calcomp page. The blade is not positioned at (0, 0) and its origin 
has no set position in relation the plotted axis. So the plot is only useful for visual ex- 
amination of the blade. 

The portion of the program which generates the Calcomp plot is coded specifically 
for the NASA Lewis system and would not work elsewhere. However, the program is 
written with the plotting code at the end. A programmer at another installation could 
easily substitute a code to obtain a plot based on the requirements of his own system. 

The coordinates for plotting are calculated and arranged in the PLOTT subroutine. 
(See COMPLETE PROGRAM LISTING.) Down to statement 310 of this routine the blade 
coordinates have been stored into two arrays: XDOWN and YACROS. The number of 
points on each blade segment have been stored into the array, NPNTS. After statement 
310, the section of code labelled PREPARE KKK AND P AND CALL CALPLT prepares 
special variables for a call on the CALPLT routine which is internal to the Lewis sys- 
tem. The subroutine CALTIT, which writes the input and output to the left of the plot, 
also uses special Lewis routines. Finally, the statement DECK CALPLT calls in the 
CALPLT routine which does the plotting. So from statement 310 of the PLOTT routine 
to the end of the coding, changes would have to be made by a programmer to get plotting 
on another system. 


Error Conditions 

Several error messages are given by the program under certain conditions. This 
section lists the error messages and explains what to do if they are encountered. 

(1) MAX THICKNESS OF SOME SEGMENT IS LESS THAN LEADING OR TRAILING 
EDGE THICKNESS OF THAT SEGMENT 

This message is printed if either of the following conditions is found on any of the blade 
segments 


TMOC < 2. 0 X RIOC 
TMOC < 2. 0 X ROOC 


The blade segments must be at least as thick as their leading- or trailing-edge thick- 
nesses. 

(2) THE SUM OF ONE PLUS THE VALUES IN THE PHOPHl ARRAY MUST BE 
GREATER THAN 0. 1 
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This message is only printed when negative input values are used in PHOPHl and the 
sum of these input values is less than -0. 9. When excessive negative cambers are used, 
the prc^ram cannot converge on its iterations to calculate individual blade segment cam- 
bers. (A negative input to PHOPHl implies that the first blade segmeht will have posi- 
tive camber, and the segment corresponding to the negative PHOPHl will have negative 
camber. This situation is permitted in the program, but is physically unrealistic. So 
the error message eliminates long iterations on bad data. ) 

(3) PROCEDURE FOR SIZING OF BLADE CAMBERS HAS NOT CONVERGED IN 
25 ITERATIONS 

The program initially calculates blade segment cambers and then corrects these cambers 
in an iteration process until an overall blade camber of DELK is obtained. Usually four 
or five iterations are required to reach a specified tolerance. The error message is 
given if convergence is not obtained in 25 iterations. This error is generally due to the 
fact that specified inputs are not geometrically compatible. This condition is most likely 
to occur when DELK is small (0° to 10°). 

In addition to the programmed error messages, computer errors (such as square 
root of negative number) are likely to occur if input values are beyond recommended 
limits. Limits within which computer errors are not likely are summarized. 


1 < N< 5 
TCHORD > 0. 
SOLID > 0. 

0. < DELK < 180. 
-90. :< KAPIN < 90. 
RADIUS > 0. 

COCl > 0. 

-1000. < PHOPHl < 1000. 
0. < GOTO < 0. 5 
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-0. 2 < LOTC < 0. 4 


0.. :< F < 10. 

0. :< TMOC < 0. 8 
0. < RIOC < 0. 4 
0. :< ROOC < 0. 4 

COMPLETE PROGRAM LISTING 

$iejoe 

HeflC CATBP 

COMMON / INPUT/N ,TCHURD , SOLID, OELK,KAFIN,PADI US, CCCK5), PI-0PH1( 5 ), 1 

1G0TC( 5 ),LOTC« 5) ,F( 5J , TMOC (5) ,RICC(5>,R0CC(S), TITLE! 12) 2 

COMMON /CUT PUT/CHORD (5 ) ,GAM(5),6AMR(5),FHIC(5),PITCH 3 

COMMON/CLPLOT/XPEN, YPEN,NX,NY, I P£N,XLAeEL{10),YLABEL(10) A 

COMMON/COMl/RK 5) ,RO( 5) ,THETA(5I ,XI(5),YI(5),XC(5),YC(5), 5 

1XCMI5),YCM15),XH5),Y1(5),X2(5),Y2(5),SLS(5),SLP(5»,XSM(5),YSM(5), 6 

1XPM( 5 ),YPM(5),KCM(5),G(5),GA(5),GAOC(5),L(5),FA(5),SINC(5), 7 

IRCI 5 ) , HC( 5 ) ,BC ( 5) ,K IC ( 5) ,KOC (5) ,PS(5),HS(5),BS(5),KIS(5),KCS(5), g 

1PHIS(5),RP!5),HP<5),BP(5I,KIP(5),K0P(5),PHIP(5),TM(5),XR(5),YR(5), -9 

1XG(5),YG(?),N0EL(5»,XX(5,1CC),YS(5,100),YP(5,100) 1C 

COMMON/COM2/GAMS,MCHORD(5) ,RSTGR(5),STGP(5),MLE(5),PTHLE(5), 11 

1THLEI5 ),MTE(5),RTHTE(5),THTEt5) ,MCL(5) ,RTHCL(5),THCL(5»,MCT(5), 12 

IRTHCTI S),THCT(5),8ETIS(5l,BeTCS(5l,BETIP(5),BETCP(5),MSPS(5,100), 13 

1RTHSPS(5,1CO»,THSPS(5,100),MSFP(5,100I ,RTHSPP ( 5 , 100 ) ,THSPP( 5, 100 ) lA 

REAL L,LOTC,NEWC,KIC,KOC,KIS,KCS,KI F,KCP,KCM, 15 

1KGS1,KC-S2,KGP1,KGP2,KGSI ,KAPIN 16 

REAL MLE,MTe,MCL,MCT,MCHCRD ,MSPS,MSPP 17 

1C CALL BLDCRO 18 

CALL IFJiNIPT 19 

CALL PLCTT 2C 

GO TO 10 21 

END 22 


SIEFTC BLCCR 

SUBROUTINE BLOCRD 1 

COMMON/INPLT/N,TCHORD ,SOLID,DELK ,KA P I N , P AC lUS , COCl ( 5 ) .PHOPHK 5 1, 2 

IG0TC( 5 ),LOTC ( 5 > ,F ( 5) , TMOC (5) ,RI CC(5) ,RC0C (5 ) ,T ITLE( 12) 3 

COMMON /OU TP UT/CHORO(5) ,GAM(5),GAMR(5),FHIC(5),PITCH A 

COMMON /COM 1/R I( 5) ,RO(5) , THETA (5) ,XI(5),YI(5),XC(5),Y0(5), 5 

1XCM{ 5 ),YCM(5),X1(5),Y1(5),X2(5I,Y2(5) ,SLS(5) , S LP (5 ) ,XSM( 5) , YSM( 5 ) , 6 

IXPM ( 5 ) ,YPM( 5) ,KCM(5) ,G(5) ,GA( 5) ,GAOC(5) ,L(5),FA(5),SINC(5), 7 
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1RC( 5),HC( 5>»BC(5) ,KIC (5) ,KCC{5) ,P5(5) ,HS(5) ,BS(5),KIS(5),K0S(5), 8 

1PH1S<5),RP15)*HP(5),BP(5) ,KIP(5) ,KCP(5)»PHIP(5),TM(5),XR(5),YR(5), S 

1XG(5I ,VG( 5)»NDEL(5) ,XX<5,I0C) ,YS(5flO0)»YP(5,lOO) 10 

HEAL L ,LOTC»NEWC,KIC,KGC»KIStKCS,KIP ,KCP,KCM,LCt 11 

IKGSlt KGS2,KGPltKGP2,KGSI »KAPIN 12 

C 13 

C READ AND PRIM INPUT 14 

C 15 

1C WRITE(6, ICCO) 16 

READ ( 5,1250 ( 11 TLE( I) ,1=1 »12) 17 

WRITE! 6,12 60 (TITLE ( I) ,1=1 ,12) 18 

READ ( 5, 1C 20) N , TCHORO , SOLI D , DE L K , K A F I N , P AD I US IS 

WR ITE( 6, 1030 ) N,TCHORD ,S0LID,DELK,KAF1N,RADIUS 2C 

READ (5,1010) (COCl (J) , J = 2 ,N) 21 

WRI TE( 6,10^0) (COCUJ) ,J=2,N) 2 2 

REAC (5,1010) (PHOPHK J) ,J=2,N) 23 

WRITE(6, 1050 ) (PHOPHU J) ,J = 2,N) 24 

REAC (5,1010 (GOTC(J) ,J = 2,N) 25 

WRilE(6,i060) (GOTO ( J) ,J = 2,N) 26 

REAC (5,1010) (LOTC (J) , J = 2 ,N) 27 

WRITE! 6, 1070) ( LO TC ( J ) , J = 2 , N ) 26 

REAC (5,1010 (F(J),J=2,N) 29 

WRITE! 6,1060) (F(J),J=2,N) 30 

REAC (5,1010) ( TMOC ( J) ,J = 1 ,N) 31 

WR I TE( 6, lose ) ( I.MOC ( J ) ,J = 1 ,N) 32 

REAC (5,1010) (RIOCU) ,J = 1 ,N) 33 

WRITE<6,1100) (RIOC ( J ) , J=1 ,N) 34 

REAC (5,1010) (ROOC(J) 35 

WRi TE ( 6, 1 110 ) (ROOC < J ) , J = 1 ,N) 36 

C 37 

C INITIAL VALUES OF CAPBER AND CHCRD 38 

C 39 

DELK = DELK/57. 295779 4C 

P ITCH= TCHORD/SOLIO A 1 

SUMC= 0. A2 

SUML= 0. 43 

SLMPHI^ 0. A4 

PH0PH1(1)= 1*0 45 

IF (N .EQ, 1) GO TO 30 46 

DO 20 J=2,N 

SUML= SUML+LOTC(J) ^8 

SUMC= SUMC+COCKJ) 49 

2C SUMPHi= SUMPHI+PHOPHl ( J) 5C 

3C factor^ 1./( 1 .-DELK^«2/24. ) 51 

SUML= EACTOR+SLML 52 

SUMC= l.+SUMC 53 

SUMPH1= 1 . + SUPPHI 54 

IF ( SUPPH I .GT*. 1) GO TO 40 55 

WR ITE ( 6, 127C ) 56 

GO TO 1C 57 

4C CF0RD(1)= TCHORD^SUMl/SLPC 56 

PHIC( J )= DELK/SUMPHI 59 

IF (N .EQ. 1 ) GC TO 60 eC 

DO 50 J=2,N 61 

PHIC(J) = PHIC( D’^'PHUPHK J) 62 

5C CH0RD(J)= CHORD (D’i'COCKJ) 63 

C 64 

C SIZING OF CTFER BLADE SEGPENT DIPENSICNS 65 

C 66 

tC ITER = 1 67 

IF (N .EQ. 1 ) GO TO 8C 68 

DO 70 J=2,N 69 
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L(JI= tcrC( J )*1CH0RD 7C 

7C G«J)= 60TC(J )*7CH0R0 71 

ec DO ?0 J=1,N 72 

TM(J)= TMDC( J)*CHORD( J) 73 

R I( J )= RIGC( J )*CHCRO( J» 7A 

RO(J) = ROOC( J)*CHORD( J) 75 

IF (2.=*RI (J).LE.TM(J).AN0.2.*R0(J».LE.TM{jn GC TO 90 7t 

RRITE(6,1120> 77 

GO 70 1C 78 

5C CONTINIE 79 

8C 

SEGMENT CENTER LINE CALCULATICNS 81 

82 

CO 100 J=lfN 82 

XKJI = RKJl 8A 

YI( J» = RKJ) 85 

XO(J) = CHORD! J)-RC(J) 86 

YO(J) = RO(J) 87 

ARG=( Y 1(J »-Y0( J 1 » /(XI ( J|-XC( Jl I 6 8 

ICC THETA(J)= -ATAN(ARG» 89 

lie DO 120 J=1»N 90 

KIC(J) = PHIC( Jl/2.-THETA( J) 91 

KOC(JI = -PHIC( Jl /2.-THETA( Jl 92 

IF (ABS(PhIC( Jl 1.LT..00C11 GC TO 120 93 

8C( Jl = (X I( J 1**2+YI ( Jl**2-X0( Jl **2-Y0( Jl *«2-2.* (XI ( J1-X0( J 1 !♦( XI 94 

1( JMY I ( J1*TAN (KIC( Jll 1 1 /2./( Y0( Jl-Yl ( J 1 + ( X I ( J 1 -X0( J 1 1 *T AN( K IC( J 1 1 1 95 

HC(J1= -XHJ l-(Yl (J1 + BC( Jl 1*TAN(KIC( Jl 1 96 

RC(J)= SORT! ( X I ( J l+HC ( Jl )*«= 2 + (YI (JM-eC( Jl 1*«21 97 

GO TO 120 98 

12C 6C(J1= 999.99999 99 

IF (PH 1C(JI.LT.C.1 BC(J1= -BC(J1 ICC 

HC(J1= -ChORD(J)/2. 101 

RC(J1= 999.99999 102 

12C CONTINLE 103 

1C4 

SEGMENT SERFAGE CALCLLATIONS 105 

1C6 

DO 27C J=lfN 1C7 

K1 = 0 1C8 

ITIR= 0 1C9 

CEL = 0.1*CHORD(J) IIC 

IF (R 1 (J 1 .EQ.O. 1 GO TO 140 111 

RURI =R0(J)/RI(J1-1. 112 

GO TO 150 112 

14C RORI= 0. 114 

15C ROMRI= R0(J1-RI(J1 115 

CRO = CHUR0(J1-2.*R0( Jl 116 

CV = ( TM( J 1-2. *R1 (Jl 1 /CHORD! Jl 117 

XCM(J1 = CH0R0(J)/2. 118 

XCMMl = XCM(J1 119 

LC= SGPT( (CH0RD(J1-RI(J1-RC(J11**2+(RI(J1-RC( J) )*»21/2.0 12C 

PC= PH1C( J 1/2.0 121 

HCM= LC*TAN(PC/2.01 122 

XC= R 1(JMLC=!'CCS(THETA(J11+HCM*S1N(THETA( Jll 123 

16C YCM(J1= R0(J1+(R1( J|-R0(J1 1*(CHCB0( J1-XCM(J1-RQ(J1 }/( CHORD! Jl 124 

l-RKJ 1-R0(J1 I 125 

IF (ABS(PH1C( Jl 1.LT..0C011 GC TO 170 126 

ALPH= 1HETA(J1+AS1N((XC-XCM( J||/IC*SIN(PC1-SIN(THETA(J1 11 127 

YCM(J1= YCM( J l+LC /COS (THETA! Jl !♦ ( SI N ( PC 1 / ( 1 . + COS ( PC 1 1 -S IN( ALPHI**2 128 

l/(( l.+tOS(ALPhl 1*£1N(PC1 1 1 129 

17C ARG= -((XCM(Jl-RI(Jll*SIN(PCl-LC*SIN(KIC(Jlll/( (YCM(J1-RI( Jll* 13C 

1SIN(PCI+LC*C0S(KIC(J1I1 131 
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KCM(J )= ATAN(ARG) 132 

C SUCTION SUKFACfc 133 

XSMIj; = XCM( J) /2.*SI MKCN(J) ) 134 

YSM(J) = VCM( J)+TM J) /2.*C0S(KCMJ) ) 135 

OS = X SMI J )*R0RI+Cf<0 136 

XMRI = XSM(J)-RI(J» 137 

XMRI2 = XSMIJ )-2.*RI (J) 138 

YMKI = VSMIJ)-RIU) 139 

YMRI2 = Y SMI J )-2.*RI IJ) 14C 

XMYM = XSM(J)**2+YSM( J) **2 141 

AAS = XSM(JM=XMRI2*YSM(J)**2*RORI**2-2.*XMRI*YMRI*YSM(J)*RORI*DS 142 

l+YSMIJ )*YMRI2*DS**2 143 

BBS = (XMYM*YMRI+RI (J)**3-3.*RI U)**2*VSM( J) ) *0S4*2-(XMYM*XMRI 144 

l+RKJ )**3-3.‘XRI (J)**2'i'XSM( J) ) * Y SM ( J ) *PORl *0S+ I XSMIJ ) f>X MR 12 *YSM(J ) 14 5 

2*RURI- XMR 1*YMRI *OS)*^I XM YM*R0P I + CHORD IJ I *C PO+RO IJ ) 4RCMRI ) 146 

CCS = (XM YM**2+RI IJ)**4-6.*fii I J) **2*XMYM) *DS*#2+XSMIJ) *XMR 12 147 

l*IXMYM*RORI+CHORD IJ)*CRO+ROIJ )*RCMRI ) **2-2 . * IXMYM+XMR I + R II J 1**3 148 

2-3 .*R I Id )**2*XSM ( J ) )*I XMYM*ROR 1 + CHORD I J ) *C PO+RO I J ) *POMR I ) * DS 149 

IF IRIIJ) .EQ.C. > GO TO 180 150 

BSIJ) = I -BBS+SQRTIBB S**2-AAS*CCS) ) / 12 .*A AS) 151 

GO TO 190 152 

18C BSIJ)= -BBS/I 2.*AAS) 153 

19C HSU) = -IXMYM*RCRI+CHCROIJ)*CRO+ROIJ)*RCMRI + 2.*YSNIJ)*RORI*BSIJ)) 154 

1/I2.*( XSMI J)*RCRI+CRQ)) 155 

SLSU) = - IXSM J) +HSI J ) ) / I VSM I J) +6SI J) ) 156 

C PRESSURE SURFACE 157 

XPMIJ) = XCM U )+TMIJ ) /2.*SI NI KCM I Jl ) 156 

YPMU) = YCMU)-TM< J) /2.*C0SIKCMIJ) ) 159 

DP = XPM( J )*R0P I+CRO 16C 

XMRI = XPM(J)-KHJ) 161 

XMRI2 = XPMIJ )-2.*RI IJ ) 162 

YMRI = YPMIJj-RIU) 163 

YMRI2 = YPMI J )-2.*Rl ( J) 164 

XMYM = XPMIJ )**2+YPMU) **2 165 

AAP = XPM(J)*XMRI2*YPMIJ)**2*RCRI**2-2.*XMRI*YMRI*YPNIJ)*RCRI*0P 166 

l+YPMI J )*YMRI2*DP**2 167 

6BP = IXMYM*YMRI+RI Ij)**3-3.*fill J)**2*YPMI J) )*DP**2-IXMYM*XMRI 166 

i+RI I J )**3-3.*Rl I J)**2*XPMI J) ) *YPMI J)*RORI*OP+ IXFMI J)*XMRI2*YPMI J ) 169 

2*ROR I-XMR 1*YMR I*DP )*I XMYM*RCR I+ChCRD I J ) *C R 0+R 0 1 J ) * R CM R I ) 17C 

CCP = IXMYM** 2+RI I J)**4-6.*RI I J) **2*XMYM) *DP**2+XPM(J)*XMR 12 171 

I* I XMYM«RORI+CHORDI J)*CRO+RGIJ)*RCMRI ) **2-2 . * I XMYM*XMR I +P II J )**3 172 

2-3 .*R 1 I J)**2*XPM I J ) )* I XMYM*RCRI+CHCRC I J ) *CP0+RC 1 J ) * PCM R I ) * CP 173 

IF IRIIJ) .EQ.C. ) GO TO 200 174 

BPU) = { -BBP-SQRT I BBP**2-AAP*CCF) ) / I2.*A AP) 175 

GO TO 210 176 

2CC BPIJ)= -BEP/I 2.*AAP ) 177 

210 FPU) = - (XMYM*RORl+CHORD I J)*CRO+ROI J) *POMRI + 2.*YPMI J) *ROR I + BPI J ) ) 176 

l/I 2.*I XPMI J )*RORI+CRO) ) 179 

SLPIJ) = -IXPMIJ)+HPIJ) )/IYPM(J)+BPIJ) ) 16C 

C CHECK FOR MAX THICKNESS POINT CCNVEFGENCE 181 

CSL = SCSI J)-SLPI J) 182 

IF ICV.EQ .C. ) GO TO 260 183 

IF lAESIUSU.LE ..CC0l*CV) GO 1C 260 184 

IF IITIR.EQ.O) GO TO 22C 185 

IF (DSL/USLM1.LT..O) Kl=l 186 

22C ITIR= ITIR+i 187 

IF IK 1 ,EQ .C) GC TO 23C 168 

GO TO 240 189 

23C XCMIJ) = XCMl J)+OSL/ABSIDSL)*OEL 19C 

GO TO 250 191 

24C XCMIJ) = XCMM1+IXCMM2-XCMM1) /Il.-CSLMl/OSU 192 

25C CSLMl = DSL 193 
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XCMM2 - XCMMl 1S4 

XCMMl = XCM(J) 195 

GO 10 160 196 

C FINAL CALCULAHONS AFTER CONVERGENCE 19? 

26C RS(J) - SQRT( (XSM( J)+HS(J) )^*2+(YSM( J)+BS( J) 198 
RP(J1 - SQRTi <XPM(J) + HP( J) )**2+ (YPM( J) + BP( J) >**2) 199 

ARG = -{R I( J )+HS( J) ) / (RI ( J)+BS( J) ) 2CC 

KIS(J ) = ATAN(ARG) 2C 1 

ARG = -(R 1(J)+FP( J) )/(RI( J)+BFC J)) 2C2 

K IP( J > = ATAN (ARG) 2C3 

ARG = -(CHORD( J)+HS(J)-RO( J) ) /(R0( J)+eS( J) ) 2CA 

KOS(J ) = ATAN( ARG) 205 

ARG = -(CHORO( J )*HP( J)-R0 ( J) ) ✓ (ROC J) +BP( J) ) 2C6 

KOP( J ) = ATAN( ARG) 2C7 

PHIS(J) = KIS(J)-KOS( J) 2C8 

27C PHlPiJl = KIP( J)-KOP( J) 2C9 

210 

LCCATION CF BLADE SEGMENTS UITH RESPECT TC CNE ANCTHER 211 

212 

GAM( 1 ) =: 0. 212 

GAMR( 1 )= 0. 21A 

IF (N ,10. 1 ) GO TO 33C 215 

DO 210 J=2,N 216 

XR(J)= CHORD( J-1)«RG( J-l)-L( J) 217 

IF (BP(J-l).LT..O) GO TO 280 218 

YR( J )= SQRT(KP ( J“1)**2-(XR( J)+HP(J-1) ) «♦£ ) -BP ( J- 1 ) - F ( J ) *G( J)-RI(J) 219 

GO TO 290 22C 

26 C YR( J ) = -S0RT(RP(J-l)**2-(XR( J)+HP(J-1) ) ^*2 ) -BP ( J-1 ) -F ( J ) ♦ G( J ) -R I ( J ) 221 

290 XG(J)= CHORD! J-l)-KO(J-i) 222 

YG(J)= -G(J) 223 

AA= 2.*C ( XG( J )-XK ( J) ) (RI ( J) +HSC J) >+ (YG( J )-YR( J ) )* ( R I( J )+BS< J ) ) ) 22A 

BB= 2.*< ( XG( J )-XR( J) ) ♦ (RI ( J) +B$( J) ) -( YG( J)-YR (J))*(RI(J)+HS(J))) 225 

CC= R 1 ( J )♦( 2.0 + RSC J)-RU J) )-(XG( J) -XR ( J) ) «*2 - ( Y G ( J ) -YR ( J ) ) ^^2 22 6 

IF (L(J).LT.O.) GO TO 300 227 

SIN GAM = (8B*CC-*AA*SQRT(AA*«2 + 8B^*2-CC* + 2) ) / ( AA*«2 + ee + *2) 226 


GO TO 310 

20C SINGAM = (BB*CC+AA*SQRT(AA^*2 + BB**2-CC**2 ))/( /iA^*2 + eB**2) 

21C GAM(J)= AR SIN( SINGAM) 

CHECK CN OVERALL ELAOE TURNING AND RESIZING CAMBERS CF BLADE SEGMENTS 
CO 320 J=2,N 

22C GAMR(J)= GAMR( J-U+GAM( J) 

23C DELKT^ K I C ( 1 ) +G AMR ( N ) - KCC ( N ) 

CIFF= OELK-DELKT 

IF (ABS(OIFF ) .LT.,001) GO TG 360 
ITER = 1 TER 4 I 
IENC= C 

IF ( M £R .GT.25 ) GO TO 350 
00 2A0 J=1,N 

2AC PHIC(J)= PHIC( J)+DIFF«PH0PH1( J) /SLMPHI 
GO TO 370 
25C HRilE(6,1120) 

GO TO 10 

RESIZING CHORDS OF BLADE SEGMENTS 

36C IEND= ] 

27C X06 = XO(N) 

YOB = >0(N) 

IF (N .EQ. 1 ) GC TO 39C 
DO 380 K=2tN 


229 

23C 

231 

232 

233 
23A 

235 

236 

237 
236 
239 
2AC 
2A1 
2A2 
2A3 
2AA 
2A5 
2A6 
2A7 
2A8 
2A9 

250 

251 

252 

253 
25A 
255 
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J= N-K-»2 256 

SING= SIN(GAM(J)) 25? 

COSC= COS(GAM{Jj) 256 

XOBB= >R{ J )+XCB*COSG+YQB*SING-RI (J) *(SING+CaSG) 259 

YOBB= YR( J )+YCB*COSG-XCa*SING+RI IJ) *(SIN6-CGSG) 26C 

XOB = >08E 261 

'8C YOB = YCBE 262 

29C IF (IENC.bg. I) GO TO AlC 263 

NEWC = EQR,T( (XOB-XKl) )**2+(yCB-YI(l) )*=»2)+RI(l)+RC(N) 26A 

CRATIO= TCHORO/NENC 265 

CO ACC J=1,N 266 

ACC CbiORD(J)= CHORD ( J )*CRATIO 267 

GO TO fO 268 

269 

OVERALL BLADE THETA AND GARMA 27C 

271 

AlC ARG= (RI( 1 )-RC(N) )/(TCHCRO-RI (l)-RO(Nn 272 

THETB = ATAN(ARG) 273 

ARG= ( yl( 1 )-YOB )/(XC8-XI (in 27A 

GAMBTE= ATAN(ARG) 275 

GAME= CAM8T8-IFETB 276 

IF (N .EU. 1) GO TO A30 277 

278 

PSEUDO- INC lOENCE ANGLES ON AFT BLADES 279 

28C 

DO A2C J=2.N 281 

ARG= (CHORD! J-1 )+HC(J-l)-RC(J-l) )/RC(J-n 282 

RH01= ^RSIN(ARG) 283 

ARG= ( ChORU( J- 1 ) +HC ( J-l)-RC( J-1 J-L ( J) )/RC ( J~l ) 28A 

RH02= ARSIN(ARG) 285 

RHO= Rb01-Rh02 286 

A2C SINC(J)= -(KIC(J)-GAM( J)-KCC( J-D-RHC) 287 

288 

ELACE SECllON CQCRCINATES AT DELX INCRERENTS 289 

C 290 

A3C DO A9C J= 1 ,N 29 1 

TEM = CHORD! J ) /2C. /lOOOC. 292 

NEXP = 0 293 

AAC NEXP = NbXP+1 29A 

TEM = lO.^TEM 295 

IF (TEM-1..LT.C.) GO TO AAO 296 

M = TEM 297 

IF (M.GE.2) GO TO A5C 298 

M = 1 299 

GO TO A70 3CC 

A5C IF (M.GE.5) GO TO A6C 3CI 

M = 2 302 

GO TO A70 3C3 

A6C M = 5 3CA 

^7C DELX = FL0AT(M)*1C.** (A-NEXP) 3C5 

NDEL(J)= CHORD! JI/DELX+1. 306 

XX! J, I ) = 0. 3C7 

NDELJ = NCEL ( J ) 3C8 

CO A9C K=l, NDELJ 309 

YS(J»K)= SOR T ! R S( J) *=^2- ( XX ! J, K) +HS ! J> ) «*2 I -eS ! J ) 310 

IF !BP!J) .LT..OI GO TO AdC 311 

yP!J,K)= SCRKRP! J)*>!‘2-!XX! J.HI+HF! jn=»«2)-EP! J) 312 

GO TO A90 313 

A8C yP! J,K )=-SQRT!RP!J)*=!'2-(XX( J,K)+HP! J) )’»«2>-8P!J) 31A 

A9C XX(J,K + 1)= XX ! J ,K ) 4-DELX 315 

C 316 

C RELATION CF SEGMENT ORIGINS TO PRINCIPAL ORIGIN 317 
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X i( 1 ) = a . 

Yl( 1) = 0, 

X 2 ( 1 ) = CHORD ( 1 ) 

Y 2 < 1 ) = C, 

IF (N .EC. ] ) GO TO 5 €G 
DO 5 CC J= 2 ,N 

TEMA = XR ( J )-R I C J )♦ (COS (GAH< J) ) + SIN (GAN ( J) n 
TEM8 = YR ( J) +R U ( SIN (GAN( J) )-CCS (GAN( J) ) ) 

TEMC = 1ENB/CGS(GAMR( J-n ) 

TEMC = lENC^SIN (GANlR( J-in 
TEME= lEMA+TEMD 
TEMF = 1 ENE*SIN(GANIR( J- 1 ) ) 

DX = T EME*COS(GAMR( J-U) 

CY = lEMF-TEMC 
Xl( J) = XKJ-U+DX 

Yl( J) = Y KJ-D-DY 

X2(J) = Xl(J)+CHORD(J)*CGS(GANR( J) ) 

5 CC Y 2 (J) = Y 1 (J)-CH 0 R 0 (J)*SIN(GANRU) ) 

CCNPLTATICN CF ACTCAL GAPS 

DO 57 C J= 2 »N 
C REAR PCR 1 ICN CF GAP 
KGSl = 100 . 

EIC ARG = - ( ( XG( J )“XR { J n^^'CCSCGAN ( J) I- ( YG (J)-YR( J) I^S IN (GAN< J) ) 

1 +R I( J ) ^HS( JM / ( { XG ( J) -XR( J) ) «SI N (GAM ( J) ) <- ( YG ( J ) - YR ( J ) )*COS(GAM( J ) ) 

2 -*RI(J HeS(J) ) 

KGS 2 = ATAN(AKG) 

IF (A 6 5 (KG 52 -KGS 1 ).LE. .on GC TC 520 
BETA = KGS 2 '-GAM(J) 

CA = ()(G( J *1 H-VG( J- 1 ) ’«^TAN(BEtA)-XR( J) ) *CCS (GAN( J )) + 

1 YR( J)«SIN(GAM( J))+Ri( J)+HS( J) 

C 0 = - (TAN (BETA )=^COS(GAM( J) J + SIN (GAM(J) ) J 
CC = (XO(J-l) + YG(J-l)*TAN(BETA)-XR(J))*SIMGAM(jn- 
1 YR( J)*CCS(GAM( J >)+RI(J) + 8 S(J) 

CD = COS( GAM( J ) )-TAN(BETA) *Si N(GAN( J) ) 

CE = CB*"^ 2 +CD ««2 
CF = 2 .^( CA*CB+CC<'CD ) 

CG = CA«^ 2 +CC** 2 -RS(J)= 4^*2 

YG(J) = (-CF + SGRT(CF4^«2“A.*CE«CG))/(2.=frCEJ 
XG( J ) = XO ( J' 1 )- ( YG ( J)-YO ( J-L > ) *TAN( BETA) 

KGSl = KGS 2 
GO TO 510 

52 C GA(J) = SORT ( ( XG( J)-XO ( J- U ) ** 2 + (YG ( . ) -YC( J -1 ) ) ’**2 ) -RC( J -1 ) 
GAOC(J)= GA( J )/CHORD( J-i) 

C FCRWARC PCRTIGN OF GAP 
KGPl = ICO. 

X GP = XR ( J ) 

YGP = S 0 R 1 (RP( J-i)** 2 -(XGP+HP(J-i ) ) **2 ) -BP ( J - 1 ) 

53 C ARG = -<XGP+-HP( J- 1 ) )/(YGP + BP( J-i )) 

KGP 2 = ATAN(ARG) 

IF (ABS(KGP 2 -KGP 1 ).LE.. 01 ) GC TC 560 
ARG = (-HS(J)-RI( J) )/(BS(J)+RI(J)) 

KGSl = ATAN(ARG)-GAM( J) 

BETA = (KGP 2 +KGSI )/ 2 . 

CL = ].4( TAN(8ETA)i=fr=i'2 

CM = 2 .=i=(BP( J -1 )-(XR{ J) + YR(J)*TAN(BETA)+HP(J- 1 ) )*TAN(BETA) ) 

CN = ( XR( J )+YR ( J) «TAN( 8 ETA)+HP( J- 1 ) ) «* 2 +BP( J-i )»* 2 -RP( J-i ) **2 

IF (BP(J- J).LT..O) GO TO 540 

YGP = (-CM + SQRT (CM** 2 - 4 .*CL«CN) ) / ( 2 . «CL ) 

GO TO 5 50 

54 C YGP = (-CM-SQPT(CM** 2 - 4 .*CL*CN) )/( 2 .«CL) 


318 
31S 
32C 
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‘50 XGP = XR( J )-( YGP-YR( J) )*TAN(BETA) 381 

KGPl = KGP2 282 

GO TO 520 282 

56C FG = SORT( (XR(J)-XGP»*^2+(YP( J)-YGP»**2> 38A 

BETA = BETA+GAP(JI 385 

CL = 1 .4( TAN(BETA)>**2 366 

CM = 2.*(BS{J)-(XI(J)+YI(J)*TAN(BETA)+HS( J»»*TAN(BETAn 387 

CN = (XI( J) + YI( J1*TAN(BETA>4H£( J)>**2 + eS( J)*<'2-RS<J)’»*2 388 

YGS = (-CM + SQRT (CM**2-4.*CL*CM ) /(2.*CL) 385 

XGS = XI( J)-(YG£-YHJn*TAN(BETA) 25C 

FG2 = SQRT((XGS-XI( J))**2+(YGS-YI(J) )**2) 251 

FG = FG-FG2 252 

57C FA(J1 = FG/GA(J» 3C2 

C 35A 

C PLT OCTPUl ANGLES IN DEGREES 355 

C 356 

*8C CELK = CELK*5T. 255775 357 

GAMB = CAM6»57. 255779 398 

TFETB = TFETB*57. 255775 39c 

00 550 J=1,N ACC 

GAM(J) = GAM(J)*57. 255775 ACl 

GAMR(J) = GAMR(J)*57. 255775 AC2 

TFETA(J) = THETA(J)*57. 255775 AC3 

SINC(J)= £INC(J)*57. 255775 ACA 

PhICCJ] = PHIC(J)*57. 255779 A05 

PFIS(J) = PHIS(J)*57. 255775 AC6 

PhIP(J» = PHIP(J>*57. 255775 AC7 

KIC(J) = KIC(J)*57. 255775 AC8 

KIS<J» = KIS( J »*57. 255775 AC5 

KIP(J) = KIP(J>*57. 255775 AlC 

KOC(J) = KOC(J)*57. 255779 All 

KOS(J» = KCS(J)*57. 295775 A12 

55C KOP(J) = KOP(J)*57. 255779 A13 

C AlA 

C CHANGE SICN CF SELECTED OLTPLIS A15 

C A16 

00 6C0 J=1,N A17 

HS( J >= -H£( J) Aie 

hC( J ) = -HC ( J » A15 

HP(J)= -HP(J) A2C 

KOS(J )= -KOS< J) A21 

KOC(J -KOC( J» A22 

KOPIJ >= -KOP( J) A23 

Y1(J)= -YKJ) A2A 

6CC Y2( J )= -Y2( J) A25 

YOB= -TOB A26 

C A27 

C PRINT OUlPLT A28 

C A25 

WRITE! 6,1CCC I A2C 

WRI TE( < ,1 lAC) N, TCHOROfPI TCH, SOLID ,OELX , KAFI N, GAMB, THeT8f X C8t YOB A21 

DO 62C J=1,N A32 

WRITE(f,l IfC ) J A23 

WRiTE((,116C) CHURO(J) ,RI <J1 ,RC( ^) , THETA! J) A3A 

WRITE!6,117C) X I ! J ) , YI ! J) , XC ! J) , YC ! J ) , XCM ! J) , YCM ! J ) A3‘ 

WRITE16,ll£0) X 1! JI ,Y1 ! J) ,X2( JJ , Y2 ! JI ,GAM! J» ,GAMR! J » A36 

WRITE!6,1150) PH I S ! J ) , RS ! J) ,H S ! J I ,BS ! J > ,K IS ! J ) , KCS ! J ) A37 

WRI7E! 6, 1200 » PH IC ! J ) ,RC ! J» ,HC ! J ) , BC ! J ) ,K I C ! J » ,KCC ! J I A38 

WR1TE!C,1210) PhIPlJ) ,RP1J) ,HP! J) ,BF!J) ,KIP!J) ,KCP! J) A35 

IF ! J .£Q. 1 ) GO TO 61C AAC 

WRITE!C, 1220 ) G ! J ) ,GA ! J I ,GACC IJ I ,L ! J ) , F ! J I ,F A i J ) ,S I NC ! J I AAl 
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tic WHI TE< t ,123C ) NDEL(J> 

NOELJ = NCEL(J) 

62C hRne(<,124C» ()(X( J,K) , YS (JtK» ,YF( JfKI , K=l.MJELJ) 

RETURN 445 

446 

FORMAT SIAlEMbNTS 447 

446 

ICCC FORMAK IHl///) 444 

ICIC FORMA 1 (EMC. 5 ) 45C 

1C2C FORMAT(IlCt5FlC.5l 451 

1C3C FORMAT(//eXtlFN,6X,6HTCHOROt6X,5FSCLIO,6X,4HOELK,6X,5hKAPINt 5X, 4 52 

16HRAU H.J/5X, I2f 5X,F8. 5 ,3X fF8.4,3XiF7. 3 f4X,F7.3 ,3Xt F8.5 ) 453 

1C4C F0KMA1(/6X,12HC/C(1) ARR A Y / ( 19X , 5 ( F 9 .5 , IX ) ) ) 454 

1C5C FORMAT(/6X,16HPHI/PHUU ARRA Y/ ( 19X ,5 (F9 . 5 ,1 X ») J 455 

1C6C FURMAT(/6X,10hG/TC ARR AY / ( 1 9X ,5 ( F9 . 5 , 1 X ) ) ) 456 

1C70 FURMAT(/6X,1CFL/TC ARR A Y / ( 1 9 X ,5 ( F9 . 5 , 1 X )) ) 457 

1C8C F0RMA1(/6X,7HF ARRA Y/ ( 1 9X , 5 ( F 9 . 5 , 1 X ) ) » 458 

1C9C FORMAT (/6X,10HTM/C ARR A Y / ( 9X , 5 ( F 9. 5 f I X ) ) ) 459 

HOC F0RMAT(/6X,1CHR1/C ARR A Y / ( 9 X , 5 ( F9. 5 . 1 X ) ) > 46C 

me FORMA T (/6X,10hRG/C ARRA Y / ( 9X , 5 ( F9. 5 , 1 X ) ) ) 461 

112C FORMAT (7//////10X,93HMAX THICKNESS CF SOME SEGMENT IS LESS THAN LE 462 

lACING CR TRAILING EDGE THICKNESS CF THAT SEGMENT) 463 

113C FORMA 1 (// /////ICX, 72HPRCC EOURE FCP SIZING CF ELACE CAMBERS FAS NOT- 464 

1 CONVERGEC IN 25 ITERATIONS) 465 

114C FORMAT ( 10Xf24HC\,ERALL BLADE P AR AMET E P S/ 14X , 1 HN , 6X, 6 FTC FORD , 7X, 466 

15HP ITCF.6X, 5HSCLID. 7X.4H0ELK, 7X,5HKAPIN,7X,4HGAMti,6X,5HTHET6, 7X, 46 7 

23FxOB,ex, 3FYOB/13Xf I2,5X,F8.5,4X,F8.5,4X,F7.4,3X,F8.3,3X,F9.4, 468 

33X, F9. «, 3X, F7. 4, 2X,F 9. 5f2X,F9. 5) 469 

115C FURMAT(///iCX,l£HBLADE SEGMENT NC. .12) 47C 

1 16C FORMAT (/ 13X, 5HCFORO .6X.2HRI , 8 > ,2 FPC , 7 X ,5H IHET A/ lOX . 4 ( F9 .5 , IX ) ) 4 7 1 

117C F0RMAT(/14X,2HX1,8X,2HYI , 8 X . 2H XC ,dX ,2HYC , 8 X , 3HXCM, 7 X, 3HYCM/ 472 

1 1CX.6 ( F9. 5,1X ) ) 473 

116C FORMA 1 ( / 14X. 2HX 1 ,8X,2hYl , 8X , 2H X2 ,8 X . 2H Y2 , 7 X » i hG AM , 7 X , 4FG AM R/ 4 74 

IlOX .61 F9. 5, IX) ) 475 

1 19C FORMA T ( / 13X.4HPHI S. 7X ,2HRS.8X,2HFS.3X,2H6S.7X,3hKIS.7X,3FKCS/ 476 

110X.6(F9.5,1X) ) 477 

12CC FORMA 1 ( / 13X. 4HPHIC .7X , 2HRC . 8 X , 2 HFC ,8 X , 2H8 C , 7X , 3 FK I C . 7X , 3HKCC/ 47 6 

llCX.6(F9.5,iX) ) 479 

12 1C FORMA I(/13X,4FPHIP,7X,2FRP,8X,2HFP,8X,2HeP,7X,3FKIP,7X,3FKCF/ 46C 

IICX. 61 F9. 5, IX ) ) 481 

12 2C FORMAT ( / 14X. IF.G, 9X, 2HGA , 7X , 4HGA CC . 7 X ,1HL,9X,1HF ,9X ,2HFA, 7X ,4HS INC/ 462 

IlOX, 7 ( F9 . 5, IX ) ) 483 

123C FORMAT (/ 14X.2HXX, 8X ,2HYS ,8X .2HYF ,5X .7HNDEL = .12) 464 

124C FORMATK 1CX.3IF9.5.1X) )) 485 

125C fORMAT(12A6) 486 

126C FORMAT ( IX , 12A6 ) 487 

127C FORMAT(///////10X.75HTHE SUM CF CNE PLUS THE VALUES IN TFE PHUPHl 488 

lARRAY MLST BE GREATER THAN C.l) 489 

END 49C 


$IEFTC IFINP 

SL8RULT1NE IFINPT 1 

COMMGN/INPLT/N.TCHORO.SCLID.OELK ,KAPIN.RAC1US.CCC1(5).PFQPH1( 5). 2 

IGOTCI 5 ).LOTC( 5) ,F (5) .TMCC(5) .RI CC(5) .ROQC { 9 ) , T I TLE ( 12 ) 3 

CCMMON/OUTPLT/CHORD(5) .GAM (5) .GAMR { 5 ) , PHI C (9 ) . P I TCH 4 

COMMON /CUM 1/R U 5 ) .RU( 5) , THETA (5 ) ,XI (5) .YI (5) .XC (5 ) .YC(5). 5 

1XCM(5),YCM(5).X1(5).Y1(5) .X2(5).Y2(5),SLS(9).SLP(5).XSM(5).YSM(5), 6 
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1XPM(5) |YP^( d) ,KCrt(5) ,G(5» ,GA(5) ,GACC (5>)tL(6),FA(5).SINC<5), 7 

1RC( 5) f FCt 5 ) f BC ( 5 ) ,KIC ( 5) tKCC (5) ,PS(5),HS(5),BS(5),KIS(5),KCS(5), 6 

iPHIS<5),RP(5)fhP{5),dP(5),KIP(5),KOP(5)fPHIP(5>,TM(5)tXR(5»tYR(5», <5 

1XG( 5 ) tVG( 5) , N0£L(5>) ,XX ( 5 , 100 » t YS (5 , 1 00 ) , YP (5 , 100 » IC 

CCMMUN/CUR2/GAf'S,MCHQRU( 5) , R S 7G fi (5 I . S T GR ( 5 ), MLE ( 5 ), RT HL E ( 5 ) , 11 

1TPLE( £ 1,HTE( 51 ,RTHTE(5) ,THTE(5) ,RCL(5) .RTHCUb ) t THCL( 5 ) t «CT( 5) , 12 

1RTHCT(£), 1hCT(£),BET1S(5),8E7CS(5) ,BEIIP(5>,BETCP(5»,MSPS(5, 100), 12 

IKThSPS (5, ICO ) , IHSPSl 5, IOC) .WSPP (5 ,100 ) ,RTHSPP ( 5 , 100 ),THSPP( 5, 100 ) 1^ 

REAL L ,L0 1C,NEV.C ,KIC ,KCC ,KI S , KCS ,KI P,KCP,KCM, 15 

1KGS1,KGS2,KGP1,KGP2,KGSI ,KAPIN IE 

REAL PL£,NfE,MCL,PiCT,MChORU,FSPS,PSPP 17 

C 18 

C 19 

C CCMPLiTATICN OF GEC')«e TRICAL INPLT FOR TANDEM BLADE, IDEAL FLOW PROGRAM 2C 

C 21 

C CFANGE SIGN CF SELECTED PARAMETEFS 22 

C 23 

CO 10 ^=1,N 2A 

KOSIJ )= -KOS( J) 25 

KOCtJ )= -KQC( J) 26 

KOP ( J )= -KQP ( J ) , 27 

Yl( J)= -YKJ) 2B 

1C Y2( J ) = -Y2 ( J ) 29 

Yaa= -YGB 2C 

C LCCATICN CF CENTERS CF LEADING EDGE CIRCLES 31 

C 32 

GAMS = (KAP IN-KICl 1 ) ) /57. 2957 79 2 2 

CO 20 0=1, N 3A 

GAMR(J) = GAMk(J)/57. 295779 35 

GAMJ = GaMS-GAMR(J) 26 

TbMl = (X1(J)-RI(1))4'C0S{GARS)-(Y1 (J)-RU1))*SIMGAMS )+PI( 1) 27 

TE.M2 = (X 1 (J)-RI ( 1) )*SIN(GAMS)+ ( Y1 ( J)-RI ( 1 ) )*CCS(GAMS ) 28 

MCL(J) = lEMl-RK J)*SIN(GAMJ)+RI (J) *CCS{GAMJ) 39 

RTFCL(J) = TEM2<-R1(J)»CCS(GARJ)+PI( J)*SIN(GAMJ) ^C 

2C THCL(J) = RTHCLUl/RADIUS 41 

C 42 

C LCCATICN CF CENTERS OF TRAILING EDGE CIRCLES 43 

C 44 

DC 20 J=1,N 45 

GAMJ = GAMS-GAMR(J) 46 

TEMl = (X2 ( J )-R I ( 1) )*COS(GAMS )-( Y2( J) -PI (1 ) ) ’('SIMGARS ) +RK 1 ) 47 

TEM2 = (X2( J )-K I ( 1) H'SIN(GAMS)<- ( Y2 ( J) -PI ( 1 ) ) *CCS (GAMS ) 48 

MCT(J) = lEMl-RCI J)*SIN(GAMJ)-RC IJ)*CCS(GAMJ) 49 

RTFCT(J) = TEM2+R0I J)*CCS(GAMJ)-RC( J)*SIN(6AMJ) 5C 

2C THCKJ) = RTriCT U) /RAOICS 51 

C 52 

C LCCATICN CF LEADING EDGES 52 

C 54 

DO 40 J=1,N 55 

MLE(J ) = MCL(J)-RKJ) 56 

RThLE(J) = RTHCL(J) 57 

4C TFLE( J ) = THCL(J) ' 56 

C 59 

C LCCATICN CF TRAILING EDGES 6C 

C 61 

CO 50 J=1,N 62 

MTE(J ) = MCT( J)+RC(J) 63 

RIhTE(„ ) = KTHCTI J) 64 

5C TFT£( J ) = THCKJ) 65 

C 66 

C LCCATICN CF LCCAL BLADE CHORDS AND STAGGERS 67 

C 6 8 
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CO 6C w=i,N 

MCHORD(J) = MTE ( J)-ML6 ( J) 

RS1GR(JI = RTHT£( J)-RTHLE(J) 

6C STGR(J) = THTt( J)-THLE(J) 

C 

C LCCAllGN CF SPLINE CLRVE ANGLES 
C 

CO 7C J^1,N 

GAMJ = GA^^S-GANR{ J) 

BETIS(J) = KiS( J)+GAMJ<^57.2<;5 779 
6ET0SUJ = KOS( J)+G AM J*57. 295779 
BETJP(j) = KIP(J)+GAMJ«57. 295779 
7C BETOP(v) = KUP(J)+GAMJ=#^57. 295779 
C 

C LOCATION CF SPLINE POINTS ON BLACES 
C 

DO SC J=1,N 

GAMJ = GAMS-GAMR(J) 

TEM 1 = (X M J )~R n 1) )=^=COS (GAMS )- (>1 U) - R I ( 1) ) *SI N (GAf'S ) I( 1 ) 

TEM2 = (Xl(J)-Rl(l))=<^SiN(GAMS)+(Yl(J)-RI(l)) «COS( GAMS ) 

NDELJ = NIDEL( J> 

CO 80 H=1,N0ELJ 

MSPS( J tK ) = TEMl + XX ( J ,KJ«=C0S(GAMJ)-YS ( J,K) *S IN( GAMJ )-MLE( J ) 
MSPP(J,K) = TEMH-XX(J,K)X^CCS(GAMJ)“YP(J,K)’{'SIN'(GAMIJ)-MLE(J) 
RTHSPS(JtK) = TEM2+XA ( J, K) ^SI N( GAMJ ) ^‘YS ( J »K ) ^S'CCS { G A M J ) - RTHL E < J ) 
RTHSPP{J,K> = TEM2+XX(J,K)*SIN(GAMJ)+YP ( J ,K) «CCS (GAMJ)-RTHLE(J) 
THSPS(J,K) = RTHSPS( J,K) /RADI IS 
EC THSPP(J,K) = RTHSPP ( J,K) /RADI IS 
C 

C PR INT CL 1PL T 
C 

VnR ITE ( e, ICCC ) 

WR I TE ( C, 1C 1C ) 

WR ITE (C,1C2C) ( JtMCHURDI J) ,STGR(.),PSTGR(J),RI(J),RC(J), 

lMLEU),THLE(JJ,KThLE(J),MIT£(J>,ThTE(J),KThTE(J),J=l,N) 

WR I Tt ( C, 1C30 ) 

WRITE(S,1CAG) (J,BET1S(JI,3ETCS(J)*BETIP(J),8ETCP(J), 

IMCL ( J ) ,THCL( J ) f KTHLL( J) tMCT ( J ) , ThCT ( J) ,KThCT ( J ) , J= 1 ,N) 

OU 9C J=1,N 

'rtR 1 TE ( 6, lose ) J 

V^R ITE ( 6 , 1060 ) 

NDLL J = NCEL ( J ) 

9C 1^RITL(6,1C?C) (MSPS(J,K),IHSPS(J5K),RTHSFS(J,K), 

1MSPP( J ,K ) , THSPP ( J,K ) ,KTHSPP ( J ,K) ,K=1 , NOEL J) 

RETLKN 

C 

C FORMAT S ] ATEMENTS 
C 

ICGO FORMAT(1H1/////10X,38HCCMPLTED INPUT FOR IDEAL FLOW P R C G RA MS/// / ) 
1010 FOKMAT(1^XtDHBLAOE, 5X ,6HMCHCRDt AX,^*HSTGR f6X,5HRSTGRf6X,2HKi,8Xt 
12 HR Of S>, 3 HMLE,6X, AH THLE,6Xt5HRTHLE»6X,3HMTE,6X, 4 hTHTE,6X,5FRTHTE) 
1020 FORMAT ( 15 X, 1 2 , 3 Xf 1 iFi C. 5) 

103 0 FORMAT ( // //14X , SHBLADE ,5X t 5HBETI S , 5 X , 5H BE T C S , 5 X » 5H B FT I P , 5X , 5FBET0P 

If 16 X f 3 HMC Lf6X,^HrHCLf6X,5HRIhCLf6Xf3FMCTf6X,4HTHCTf6X,5FRTFCT) 

104 0 FORMAT U5XfI2,3Xf4F10.5fi0X,6F10.5) 
lOfO FORMA T (// y iCX , 1 6HBLADE SEGMENT NO. fI2) 

1060 FORMAT (/14X,4HMSPSf 6X,5HTHSPSf4X,6HRTHSPS,15X,4HMSPP,6X , 

J5FThSPPf AXfEHRTHSPP) 

1070 FORMA l((lCX,3FlC.5flCXf3FlC.5)) 

END 


69 

7C 

71 

72 

73 

74 

75 

76 

77 
76 
79 
8C 
81 
82 
63 
84 
65 
86 
67 
88 
69 
9C 

91 

92 

93 

94 

95 

96 

97 
9 6 
99 

ICC 

ICl 

1C2 

103 

104 
1C5 
106 
1C7 
108 
1C9 
lie 
1 1 1 
112 

113 

114 

115 

116 

117 

118 
119 
12C 
121 
122 

123 

124 

125 

126 

127 

128 
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$IEF1C PLT 


SLBROUINE PLCIT 1 

CCMPON /INPLT/Nf TChOROf SOLlOfCELK .KAPIN .RADIUS t CCCl ( 5 ) , PHOPH 1( 5 ) , 2 

IGOTCI 5 J,LOTC( 5),FI 5) ,TMCC (5) ,RI CC(5) ,RCCC (5 ) .TITLE! 12) 3 

COMMON /GUIPLT/CHORD (5 » .G AM ( 5 ) ,G AFR (5 ) . FHIC ( 5 ) , P IT CH A 

CCMMON/CLPLOT/XPEN.YPEN.NX.NY.IPEN.XLAeELdO) ,Y LABEL! 10) 5 

COMMON /COMl/R I ! 5 ) .RO!5 ) .THETA !5 ) ,X I I 5 ) .YI 15) . XC ! 5 ) . YC ! 5 ) , 6 

IXCM! 5) ,YCM!5I .X 1! 5) .Y1 !5) .X2!5) ,Y2I5) .SLS!5) ,SLP!5 ) .XSMI5) ,YSM! 5). 7 

IXPM! 5 ) ,YPM!5) ,KCM!5) .G!5) .GA! 5) ,GAOC !5).L!5).FA!5),SINC!5), 8 

IRC! 5 ) , K! 5) .bC! 5) .KIC !5) ,KCC!5) ,FS !5) .HS !5) .BS !5 ) . K IS ! 5 ) , KCS ! 5 ) , 5 

lPHIS!5l,RP!5).hP!5).BP!5).KIP!5).K0P!5),PHrP!5).rM!5).XR!5).YR!5), 10 

lXG!5),YGlf).NUEL!5).XX!5»lCC),YS!5,lOO) ,YP !5 .100) 1 1 

C0MMJN/C0M2/GAMS.MCH0RD!5) ,PSTGR!5).STGR15).MLE!5).PTHLE!5), 12 

ITHLE! 5 ).M7E! 5) ,RIHI£!5) .THT£!5) .MCL !5) .R THCL ! 5 ) . THCL ! 5 ) . MCT ! 5 ) . 13 

IRTHCT! 5).THCT!5),BETIS!5).bETCS!3),8ETIP!5).BETCP!5).MSPS!5,lOO). lA 

IRThSPS !£, ICO) ,1HSPS!5. 100) .MSFP!5,100) .RTHSPP ! 5 , 100 ) .ThSPP ! 5 , 100 ) 15 

DIMENS ION X1S15).Y1S!5),X1P!5),Y1P!5) .X2S!5),Y2S!5).X2P!5),Y2P!5), 16 

INDELS! 5 ),NDELP! 5) ,M! 5) ,N2 !5) .NPMS15) ,XCRX!5) .YCRY!5 ) .XIX ! 5). 17 

lY lY! 5 ) .X0X15) .YCY15) .XSI5.1CC) .XP(5 .100) .XSX15 ,100 ) .YSY15, 100), IE 

IXPX! 5, ICC ) ,YP Y! 5, 100) ,XI XC!5,10C) ,YI YC15 ,100) ,XCXC !5, 100), 19 

lYOYC! 5 ,10C ) , XULViN! 20C0) , YACRC S !2 COC ), X 1 EMP ! 1000 ) ,YT EMP ! 1000 ) , 2C 

IKKK! 25 ),P 125) 21 

EQUIVALENCE 1 X S ! 1 , 1 ) , M SP S ! 1 , 1 ) ) , ! XP ! 1 , 1 ) , R THSPS ! 1 . 1 ) ) , 22 

llXSX!i,l),THSPSll,l)),!YSY!l,l),MSPP!l.l)), 23 

IIXPX! 1 , 1) ,R IHSPP! 1,1) ) , !YPY !1 ,1 ) ,ThSFF !1 ,1 ) ) 2 A 

REAL L ,L0 IC,NE)>C ,KIC ,KCC ,K1 S, KCS ,KI F ,KCP,KCM, 25 

1KGS1,KGS2,KGP1,KGP2,KGSI ,KAP1N 26 

REAL MLE,MTE,MCL,MCT,MCHCRO ,MSFS ,MSPF 2? 

28 

29 

CALCULATICN CF INPUT FOR CALCCMP PLCTTEP 3C 

31 

32 

PLT angles IN RACIANS 33 

3A 

PI = 3.1A159265 35 

DO 10 0 = 1, N 36 

KISIJ) = K IS! J ) /57. 295779 37 

KIPIJ) = KIP!J)/57. 295779 36 

KOSIJ) = KOS(J)/57. 295779 39 

1C KOPIJ) = KCPlJ>/57. 295779 AC 

A1 

OVERALL BLADE SIZE A2 

A3 

YMAX = RS! 1 )-BS 1 1 ) AA 

YMIN = YOe-RO(N) A5 

XMAX = X08+R01N) A6 

XMIN = 0. A7 

DX = XMAX-XMIN AE 

OY = YMAX-YMIN A9 

5C 

LCCATIGN CF OVERALL BLADE ORIGIN K I 1H RESPECT TO FLCT CPIGIN 51 

52 

XT = CX/18. '2 

YT = YMAX+OX/9. 5A 

c c 

LOCATION OF POINTS WHERE LEADING ANC TRAILING EDGE RADII MEET 56 

BLADE SURFACES 

5E 

DO 2C J=1,N 59 

XISIJ) = XIIJ )-RI !J)*SIMKISI J) ) 6C 
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YlSiJJ = YI( J )*RI ( J»*CCS(KI S( J) > 61 

X1P(JI = Xl( J)*RI(J)*SIN(KIP(J> ) 62 

YIP(J) = YI<J»-RICJ)*CGS(KIP( J) ) 63 

X2S(J ) = CHORD ( J»-RC( J» *( I .+SIN <KCS ( J) ) > 6^ 

Y2S(J) = YO(J>+RO(J»*CCS(KCS( J» » 65 

X2P(J) = CHORD(J»-RO(J»*(l.-SIK(KCP(J) )» 6C 

2C Y2P(J) = YO( J)-RC(J)*COSiKCP( J) ) 61 

C 68 

C ELIMINATICN OF XX.YSt AND YP PCIMS ^CT CN THE BLADE SURFACES CS 

C 70 

DO nc J=1,N 71 

NOELJ = NDEL(J) 72 

C SLCllQN SERFAGE 72 

M = 1 74 

CO 20 K=1,NUELJ 75 

IF (XX (J,K I.GT.XISC J) 1 GO TO 40 76 

3C CONTINLE 77 

4C XS( J, 1 ) = XlSU ) 78 

YS( J, 1 ) = Y1S( J) 79 

KK = K fit 

DO 50 1<=KK, NOELJ 

IF (XX (J,K ).GT.X2S( J) ) GO TC 60 62 

M = M+ 1 87 

XS( J,P ) = XX( J,K) 8^ 

5C YS« J,M ) = YS( J,K) 85 

6C M = ] 86 

XS(J,H) = X2S(J> 81 

YS( J,P ) = Y2S( J) 86 

NDELS( J ) = M 89 

C PRESSURE SURFACE 

M = 1 91 

CO 7C K=1,NDELJ 92 

IF (XX (J,K ).GT.X1P( J) ) GO TC 80 92 

1C CONTINLE 94 

€C XP( J, 1 ) = XIP( J» 95 

YP( J, I 1 = YiP( J) 96 

KK = K 97 

DO 90 K=KK,NDELJ 98 

IF (XX (J, K 1 .GT.X2PI J) ) GO TC 100 99 

M = M* 1 ICO 

XP( J,P 1 = XX(J,K) 101 

9C YP( J,M ) = YP ( J.K) 102 

ICC M = M+ 3 103 

XP( J,M ) = X2P(J» lOA 

YP( J,M ) = Y2P( JJ 105 

lie NDELP(J) = M 106 

C 107 

C LCCATICN CF LOCAL BLADE ORIGINS l«I TH RESPECT TC PLOT CRIGIN 1C8 

C 109 

CO 12C J=],N 110 

XORX( J)=X1(J)+XT 111 

12C YORY(J) = YT-YKJ3 112 

C 117 

C LOCATION CF ELADE SURFACE CCCRDI NATES KITH RESPECT TC FLCT ORIGIN 114 

C 115 

DO 15C J=1,N 116 

SING = SIMGARRIJl) 117 

COSG = C0S(GAMR(J)1 118 

NDELJ = NOELSIJI H9 

DO 13C K=l, NOELJ 120 

XSX(JtH) = XORXI JH-XS(J»K)*CCSG+YS(JtK»*SING 121 
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13C YSY(J,H) = YGRY( J)+XSC JfK)*SING-YS(J,K)*CGSG 122 

NDELJ = MCELP(J) 123 

DO lAC K=l, NOELJ 12^ 

XPX(J,K) = XORX(J) + XP( J,K)*CG5G+^iP(J,K)»SING 125 

14C YPY(JfK) = YORY(J)+XP ( J,K)#SING->[P( JtK)*COSG 126 

15C CONIIME 127 

128 

LCCATION CF LEADING AND TRAILING EDGE CIRCLE CENTERS WITH RESPECT TO 129 

PLCT ORIGIN 130 

131 

DO 16C J=1»N 132 

SING a MMGAMK(J) ) 133 

COSG = COS(GAMR(jn 13A 

XIXIJ) = XGRX (J )+RI ( J)«( SING + CGSG) 135 

YIY(J) VCKY( J)<-RI ( J)*(S1NG-CCSG> 136 

XOX( J ) = XQRX( J)+CHORD(J)*CC SG + RC(J)^< SING-COS G) 137 

16C YUY(J) = ^fQRY(J)+CHORO(J)«SING-RC( J)«(SING + CCSG) 138 

139 

LCCATION CF ELADE SURFACE POINTS ARCLNC LEADING EDGE l^C 

lAl 

CO 17C U=liN 1A2 

ANGi = PI-KISI J )+KlP( J ) 1A3 

NKJ) = ANGl/.l lAA 

ANGI = PI/^.-GAMR( J)+KIS( J) 1A5 

NIJ = NKJ) 1A6 

CC17CH=:1,NIJ lAT 

ANG 1 = ANG1+. 1 1A8 

XIXC(J,K) = XIX( J)+RI ( J)*CCS(ANGl) 1A9 

17C YIYC(J,K) = YIY( J)-RI ( J)*SIN( ANGI) 150 

151 

LCCATION CF BLADE SURFACE POINTS APCCNO TRAILING EDGE 152 

153 

CO 18C J=1,N 15A 

ANG2 = PMKCSC J)-KUP(J) 155 

N2( J ) = ANG2/.2 156 

ANG2 = Pl/2.-GAMR( J)+KCP( J) 157 

N2J ^ N2( J ) 15E 

DO 18C K=1,N2J 159 

ANG2 = ANG2+.2 16C 

XOXC(J,K) = XCX( J)-RC( J)*CCS(ANG2) 161 

18C YQYC(JfK) = YGYU)+RC( J)=<'SIN( ANG2) 162 

163 

STORE ELACE SLRFACE POINTS INTC >0CV<N AND YACRCS 16^ 

165 

M = C 166 

DC 23G J= IfN 167 

NPNTS(v) = NDELSI J)+NDELP( J)+M ( J) + N2(J) 168 

C PRESSURE 5LREACE 169 

NDELJ = NCELPl J ) 17C 

DO 19C K=l, NDELJ 171 

M = M+1 172 

XDO^NIM H YPV(J,K) 172 

19C YACROS(M) = XP>(J,K) 174 

C TRAILING ECGE 175 

N2J = N2 ( J ) 176 

DC 2CC K=1,N2J 177 

M = M+ ) 178 

XOORN(R) = YOYC(JrK) 179 

2CC YACROS(R) = XCXC(J»K) 18C 

C SUCTION SURFACE 181 

NDELJ = NDELS(J) 182 

CO 21C K=l, NDELJ 183 
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MM = NCfcLJ-K+1 le^ 

M = M + ] 185 

XDQk(N(M) = YSY(JfMM) 186 

21C YACROS(MJ = XSMJ.MM) 167 

C LfcACINC ECGC 186 

MJ = M( J ) 169 

CO 22C K=1,N1J 19C 

M = M + ] I 9 1 

XOUk>N(N) = YIYC(J,K) 192 

22C YACRUS(N) = XIXClJ,K) 1<;2 

23C CCNTINLE 19A 

195 

RCTATt BLAOES TO NORMAL CASCADE SE111NG ISt 

197 

DO 2AC 1=1, M 198 

XTEMB(l) = ( YACPOSl I )-Xl X ( 1 ) ) =»CCS(GAMS> + ( XDCWM I )-Y IY< 1 ) ) 199 

1*SIN( GAMS I 2CO 

2^C YTEMP(l) = - ( YACRLSl I )-Xl X (1 ) >*SI MGAMS1 + ( XOCKM n-YIY ( 1 ) ) 2C1 

1*CUS(GAMS) 202 

2C3 

FINE MAXIMUM AND MINIMUM LIMITS OF PLOT, AND SHIFT BLADES 2CA 

20 5 

XMIN= C. 2C6 

CU 25C 1=1, M 2C7 

25C XMIN= AM I A 1( XM IN , X TEMP ( I ) ) 2C6 

XMAX= C. 2C9 

CC26CI=1,M 21C 

26C XMAX= AMA X 1 ( XMAX , X TE MP ( I n 211 

YMilN = C. 2 12 

C027CI=1,M 213 

27C YMIN = AM INK YM IN ,YTE MP< I ) ) 21'^ 

YMAX = C. 215 

DO 28C 1=1, M 216 

28C YMAX = AMAX1(YMAX,YT£MP(I ) ) 217 

0X= XMAX-XM'IN 218 

DY= YMAX- YMIN 219 

XT= -XM IN +CX/ 18. 22C 

YI=-YMIN+CX/9. 221 

CO 29C 1=1, M 222 

XCLmNM)=YlEMP(I) + YT 22’ 

29C YAUROS(I) = XIEMiP(I) + XT 22^ 

225 

CLPL1CA1E BLADES FOK CASCADE EFFECT 226 

227 

M M = M 2 P 8 

CO 3CC K=1,MM 229 

Ml = MM-tK 23C 

XDCkN(M) = XDCPN ( K ) +PI TCH 2’1 

2CC YACKOS(M) = YACROS(K) 232 

DO 31C J=1,N 233 

31C GAMK(J) = GAMK ( Jl+57.295 779 23A 

235 

PREPARE KXK AND P AND CALL CALPLT 236 

237 

KKK( 1 ) = A 238 

KKK ( 2 ) = C 239 

KKK( 2 ) = 2*N 2AC 

KKK ( A ) = 1 241 

CO 22C J=1,N 242 

22C KKK(J-*5) = NPNTS(J) 24’ 

DO 22C J=1,N 244 

K = J+5+N 245 
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KKKiK ) 

= NPNTS(J) 

2A6 

P ( 1 ) = 

3.C 

2^7 

P(2) = 

9, /CX*(OY+P ITCH ) + 2. 

2^6 

P( 3 ) = 

C . t 

2A9 

P( ^ ) = 

DY-»PITCH+2./S.*0X 

25C 

P( 5 i = 

10 .C 

25 1 

P(6 ) = 

C.C 

25 2 

P( 7 ) = 

10 

2 5 3 

P{ 8 ) ^ 

10 .C 

25^ 

P( 9 ) = 

0. 

255 

P( 1C) = 

c . 

256 

P( 1 1) = 

0 . 

2 5 7 

P( 12 ) = 

0 . 

258 

P( 13 ) = 

c . 

259 

P(l^) = 

90. 

26C 

NX = -] 


261 

NY = ] 


2 6 2 

DATA XLABEL!l)/lri / 

26 3 

DATA YIABEL(1)/1M / 

264 

CALL CALPLT! XDCNN ,YACRCS ,KKK, F ) 

265 

RETLRN 


266 

END 


267 


$ItnC CT ITLE 


SCSRCL 1 INE CAL 1 1 T 

CUMf-QN /INFLI/N.TCHURO , SCL 10 , D E L K , KA P I N , PA D lUS 
1G0TC( 5 ),L0TC15),F(5), TMCC (5) ,RICC(5) ,PCOC (5) t 
COME'CN/OOTPCT/CHORO(5> ,GAP(5) ,G A PR ( 5 ) , PHI C ( 5 ) 
COM PUN /CL PLIJT/>P£NfYPEN»NX,NY»IPENt)<LAeEL(10) 


,CCCI(5),PE0PH1(5), 
TITLE( 12) 

, PITCH 
,Y LABEL! 10 ) 
,TITL5(5) ,TITL6(7) 


ItHTC 


Lf fH/TC 


3) /6HT0TAL 

,6h 

TC ,3HTAL/ 



7 

5) /6riCHCK0 

t6H 

CA ,6HPBER 

,6H 

SCL,5HICITY/ 

8 

5) /6HPI TCH 

,6H 

KA ,6HP IN 

♦ 6H 

RAD, 3H lUS/ 

9 

7) /6HC/CI 1) 

,6H 

PHI ,6H/PhI !l 

>6 H ) 

G/, 

1C 

2H F/ 





1 1 

5) / 6HTP./C 

♦ 6H 

P ,6H l/C 

,6H 

RC,2H/C/ 

12 

7)/6H C 

,6H 

F ,6HHI 

,6H 

GA, 

13 


16HM f EH GAPR/ 

CALL S YMBCL ( - t .C,S.5»C.C8fTIUE tG.0t72) 
CAL L S YMdCL 1-6.C,8.75,C.15,TITL1,0.0,15) 
Call SYMBCL(-£-C,8.5»C.15fTITL2,C.C»29) 
CALL SYMBOL!- 6 .C , 7.AfC.15fTITL3,C.0t27) 
CALL S YNBLL l-£.C,t.3»0.15»TnL''+,C.0,3a) 
CALL SYMBOL (-6 . C,A.6,C.15,TITL5 ,C.0f26) 
CALL SYNBCL!-e.C,2.7fC.15,T[TL6,C.Cf36) 
CALL NL)MBEk|-6 .C , 8. 3, C.12 .TCHCRO ,0.0,5) 
Call n umb eri — a.t , 8*3, c. 12 ,oelp,c*o,a) 
CALL NUPBERl-3.1,8.3,C.li,SCLID,C.O,A) 
CALL NUMBEK|-6.C,?.2,C.12,PITCH,C.C,5) 
CALL NUPBER1-A.7,7.2,C.12,KAPIN,C.0,A) 
CALL N UP»B ER!” 3.1,7.2,C.12,PA01LS,0*0,5) 
IF !N .EC. 1 ) GO TO 2C 
YYY= e .3 
DO K ^=2,N 
YYY = YYY-C.2 

CALL NUMBER!-6.C,YYY,C.12,CCCl!vi),0.0,A) 
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CALL NUMBER(-A. e,YYY,C. 12, PHCPHKJ) ,0.0,4) 33 

CALL lYl>ldEft(-'.2,YYY,C.12,GCTC(J),0.0,4) 34 

CALL NUMi3ER(-Z.3,YYY,C.12,LCTC( J),0.0,4) 35 

IC CALL NUMbER(- 1 .5, YYY,C.12 ,e ( J » ,0.0,4) 36 

ZC YYY= 4.6 27 

CO 2C J=1,N 28 

YYY = YYY-C.2 39 

CALL NUMBER<-6.C,YYY,0.12,1PCC< J) ,0.0,4) 4C 

CALL NUMbER(-4.e,YYY,C.12,RlCC(J),0.0,4) 41 

3C CALL NUM8tR(-3.2,YYY,C.12,RCCC(J),C.0,4) 42 

YYY= 2.7 43 

CO 4C J=1,N 44 

YYY = YYY-C.2 45 

CALL NUMBER(-6. C, YYY, C. 12, CHCPD(J) ,0.0,5) 46 

CALL NUMBERC-4. 7,YYY,C.l2 ,PHIC< J),0.0,4) 47 

CALL NUM8ER(-2.2,YYY,C.12,GAP(J) ,0.0,4) 48 

4C call Nl)M3ER<-1.9,YYY,0.12,GAPR(J) ,0.0,4) 49 

RETLRN 5C 

END 5 1 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, June 22, 1970 
126-15. 
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